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We study the processes e + e~ — » K + K~ 7r + 7r~ 7, K + K ~7r°7r°7 and K + K~ K + K~ 7, where the 
photon is radiated from the initial state. About 34600, 4400 and 2300 fully reconstructed events, 
respectively, are selected from 232 fb _1 of BABAR data. The invariant mass of the hadronic final state 
defines the effective e + e~ center-of-mass energy, so that the K + K ~7r + 7r~7 data can be compared 
with direct measurements of the e+e~ -> K + K~ tt + tt reaction; no direct measurements exist for 
the e + e~ — > K + K~tt n or e + e~ — > K + K~ K + K~ reactions. Studying the structure of these 
events, we find contributions from a number of intermediate states, and we extract their cross 
sections where possible. In particular, we isolate the contribution from e + e~ — > <^>(f020)/o(980) 
and study its structure near threshold. In the charmonium region, we observe the J/tp in all three 
final states and several intermediate states, as well as the tp(2S) in some modes, and measure the 
corresponding branching fractions. We see no signal for the Y(4260) and obtain an upper limit of 
^(4260)^+.- ■ r^e < 0.4 eV at 90% C.L. 
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PACS numbers: 13.66.Bc, 14.40.Cs, 13.25.Gv, 13.25.Jx, 13.20.Jf 



I. INTRODUCTION 

Electron-positron annihilation at fixed center-of-mass 
(cm.) energies has long been a mainstay of research in 
elementary particle physics. The idea of utilizing initial- 
state radiation (ISR) to explore e + e~ reactions below 
the nominal cm. energies was outlined in Ref. [1], and 
discussed in the context of high-luminosity cf) and B fac- 
tories in Refs. [2-4]. At high energies, e + e~ annihilation 
is dominated by quark-level processes producing two or 
more hadronic jets. However, low-multiplicity exclusive 
processes dominate at energies below about 2 GeV, and 
the region near charm threshold, 3.0-4.5 GeV, features 
a number of resonances [5]. These allow us to probe a 
wealth of physics parameters, including cross sections, 
spectroscopy and form factors. 

Of particular current interest are the recently observed 
states in the charmonium region, such as the F(4260) [6], 
and a possible discrepancy between the measured value 
of the anomalous magnetic moment of the muon, — 2, 
and that predicted by the Standard Model [7] . Charmo- 
nium and other states with J PC = 1 can be observed 
as resonances in the cross section, and intermediate states 
may be present in the hadronic system. Measurements 
of the decay modes and their branching fractions are im- 
portant in understanding the nature of these states. For 
example, the glue-ball model [8] predicts a large branch- 
ing fraction for y(4260) into 07T7T. The prediction for 
g^ — 2 is based on hadronic-loop corrections measured 
from low-energy e+e~ -^hadrons data, and these dom- 
inate the uncertainty on the prediction. Improving this 
prediction requires not only more precise measurements, 
but also measurements over the entire energy range and 
inclusion of all the important subprocesses in order to un- 
derstand possible acceptance effects. ISR events at B fac- 
tories provide independent and contiguous measurements 
of hadronic cross sections from the production threshold 
to about 5 GeV. 

The cross section for the radiation of a photon of 
energy E 1 followed by the production of a particular 
hadronic final state / is related to the corresponding di- 
rect e + e~ — > / cross section cr/(s) by 



da 7 f(s, x) 
dx 



= W(s,x)-a f (s(l-x)) , 



(1) 



2E-. 



where y/s is the initial e + e cm. energy, 
is the fractional energy of the ISR photon and E c 
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\J s(l — x) is the effective cm. energy at which the fi- 
nal state / is produced. The probability density function 
W(s, x) for ISR photon emission has been calculated with 
better than 1% precision (see e.g. Ref. [4]). It falls rapidly 
as E 1 increases from zero, but has a long tail, which com- 
bines with the increasing <7/(s(l — x)) to produce a sizable 
cross section at very low E c . m .. The angular distribution 
of the ISR photon peaks along the beam directions, but 
10-15% [4] of the photons are within a typical detector 
acceptance. 



Experimentally, the measured invariant mass of the 
hadronic final state defines E c m . An important feature 
of ISR data is that a wide range of energies is scanned 
simultaneously in one experiment, so that no structure 
is missed and the relative normalization uncertainties in 
data from different experiments or accelerator parame- 
ters are avoided. Furthermore, for large values of x the 
hadronic system is collimated, reducing acceptance issues 
and allowing measurements at energies down to produc- 
tion threshold. The mass resolution is not as good as a 
typical beam energy spread used in direct measurements, 
but the resolution and absolute energy scale can be mon- 
itored by the width and mass of well known resonances, 



such as the J/ip produced in the reaction e" 1 



J/V>7- 



Backgrounds from e + e~ — > hadrons events at the nomi- 
nal y^I and from other ISR processes can be suppressed 
by a combination of particle identification and kinematic 
fitting techniques. Studies of e+e~ — > and sev- 

eral multi-hadron ISR processes using BABAR data have 
been reported [9-12], demonstrating the viability of such 
measurements. 



The K + K 7r + 7r final state has been measured di- 
rectly by the DM1 collaboration [13] for ^/s < 2.2 GeV, 
and we have previously published ISR measurements of 
the K+K-tt+tt- and K + K~K + K~ final states [11] for 
^c.m. < 4.5 GeV. We recently reported [14] an updated 
measurement of the K + K~it + it~ final state with a larger 
data sample, along with the first measurement of the 
K + K ~-K°n° final state, in which we observed a struc- 
ture near threshold in the <f>fo intermediate state. In 
this paper we present a more detailed study of these two 
final states along with an updated measurement of the 
K + K~ K + K~ final state. In all cases we require detec- 
tion of the ISR photon and perform a set of kinematic 
fits. We are able to suppress backgrounds sufficiently to 
study these final states from their respective production 
thresholds up to 5 GeV. In addition to measuring the 
overall cross sections, we study the internal structure of 
the events and measure cross sections for a number of 
intermediate states. We study the charmonium region, 
measure several J /ip and ip(2S) branching fractions, and 
set limits on other states. 



7 



II. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric energy e + e~ 
storage rings. The total integrated luminosity used is 
232 fb"\ which includes 211 fbT 1 collected at the r(4S) 
peak, \fs = 10.58 GeV, and 21 fb -1 collected below the 
resonance, at y/s = 10.54 GeV. 

The BABAR detector is described elsewhere [15]. Here 
we use charged particles reconstructed in the track- 
ing system, which comprises the five-layer silicon vertex 
tracker (SVT) and the 40-layer drift chamber (DCH) in 
a 1.5 T axial magnetic field. Separation of charged pi- 
ons, kaons and protons uses a combination of Cherenkov 
angles measured in the detector of internally reflected 
Cherenkov light (DIRC) and specific ionization measured 
in the SVT and DCH. For the present study we use a kaon 
identification algorithm that provides 90-95% efficiency, 
depending on momentum, and pion and proton rejection 
factors in the 20-100 range. Photon and electron energies 
are measured in the CsI(Tl) electromagnetic calorimeter 
(EMC). We use muon identification provided by the in- 
strumented flux return (IFR) to select the final 
state. 

To study the detector acceptance and efficiency, we 
use a simulation package developed for radiative pro- 
cesses. The simulation of hadronic final states, includ- 
ing K + K-Tr + n--f, K + K-TT°n°-f and K + K~ K + K~^, 
is based on the approach suggested by Czyz and 
Kuhn[16]. Multiple soft-photon emission from the initial- 
state charged particles is implemented with a structure- 
function technique [17, 18], and photon radiation from 
the final-state particles is simulated by the PHOTOS 
package [19]. The accuracy of the radiative corrections 
is about 1%. 

We simulate the K + K~irir final states both accord- 
ing to phase space and with models that include the 
0(1020) -> K+K~ and/or / (980) -> tttt channels, and 
the K + K~ K + K~ final state both according to phase 
space and including the <f> — > K + K~ channel. The 
generated events go through a detailed detector simula- 
tion [20] , and we reconstruct them with the same software 
chain as the experimental data. Variations in detector 
and background conditions are taken into account. 

We also generate a large number of background pro- 
cesses, including the ISR channels e + e~—t 7r + 7r~7r + 7r~7 
and 7r + 7r~7r°7r°7, which can contribute due to particle 
misidentification, and ^777, 4>ir°-f, n + ir~ir°j, which have 
larger cross sections and can contribute via missing or 
spurious tracks or photons. In addition, we study the 
non-ISR backgrounds e + e~^qq (q — u, d, s, c) generated 
by JETSET [21] and e+e" t+t" by KORALB [22]. 
The contribution from the T(45) decays is found to be 
negligible. The cross sections for these processes are 
known with about 10% accuracy or better, which is suf- 
ficient for these measurements. 



III. EVENT SELECTION AND KINEMATIC FIT 

In the initial selection of candidate events, we con- 
sider photon candidates in the EMC with energy above 
0.03 GeV and charged tracks reconstructed in the DCH 
or SVT or both that extrapolate within 0.25 cm of the 
beam axis in the transverse plane and within 3 cm of 
the nominal collision point along the axis. These criteria 
are looser than in our previous analysis [11], and have 
been chosen to maximize efficiency. We require a high- 
energy photon in the event with an energy in the initial 
e+e~ cm. frame of E 1 > 3 GeV, and either exactly four 
charged tracks with zero net charge and total momen- 
tum roughly opposite to the photon direction, or exactly 
two oppositely charged tracks that combine with a set 
of other photons to roughly balance the highest-energy 
photon momentum. We fit a vertex to the set of charged 
tracks and use it as the point of origin to calculate the 
photon direction. Most events contain additional soft 
photons due to machine background or interactions in 
the detector material. 

We subject each of these candidate events to a set of 
constrained kinematic fits, and use the fit results, along 
with charged-particle identification, both to select the fi- 
nal states of interest and to measure backgrounds from 
other processes. We assume the photon with the highest 
E 1 in the cm. frame is the ISR photon, and the kine- 
matic fits use its direction along with the four-momenta 
and covariance matrices of the initial e + e~ and the set 
of selected tracks and photons. Because of excellent res- 
olution for the momenta in the DCH and good angular 
resolution for the photons in the EMC, the ISR photon 
energy is determined with better resolution through four- 
momentum conservation than through measurement in 
the EMC. Therefore we do not use its measured energy 
in the fits, eliminating the systematic uncertainty due to 
the EMC calibration for high energy photons. The fitted 
three-momenta for each charged track and photon arc 
used in further kinematical calculations. 

For the four-track candidates, the fits have three con- 
straints (3C). We first fit to the 7r + 7r~7r + 7r~ hypothesis, 
obtaining a xlw If the four tracks include one identified 
K + and one K~, we fit to the K + K~tt + tt~ hypothesis 
and retain the event as a K + K ~ir + n~ candidate. For 
events with one identified kaon, we perform fits with each 
of the two oppositely charged tracks given the kaon hy- 
pothesis, and the combination with the lower Xkktt+tt- 
is retained if it is lower than xln ■ If the event contains 
three or four identified K^, we fit to the K + K~ K + K~ 
hypothesis and retain the event as a K + K~K + K~ can- 
didate. 

For the events with two charged tracks and five or more 
photon candidates, we require both tracks to be identified 
as kaons to suppress background from ISR 7r + 7r~7r°7r° 
and K ± K®tt t events. We then pair all non-ISR pho- 
ton candidates and consider combinations with invariant 
mass within ±30 MeV/c 2 of the ir° mass as ir° candi- 
dates. We perform a six-constraint (6C) fit to each set of 
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20 40 . 60 

X(2K2Jt) 

FIG. 1: Distribution of x from the three-constraint fit for 
K + K~ tt + tv~ candidates in the data (points). The open his- 
togram is the distribution for simulated signal events, normal- 
ized as described in the text. The cross-hatched (hatched) his- 
togram represents the background from non-ISR events (plus 
that from ISR 4ir events), estimated as described in the text. 

two non-overlapping 7r° candidates plus the ISR photon 
direction, the two tracks and the beam particles. Both 
7T° candidates are constrained to the 7r° mass, and we 
retain the combination with the lowest Xkk • 

IV. THE K+K--K+-K- FINAL STATE 
A. Final Selection and Backgrounds 

The experimental Xkktt+tt- distribution for the 
K + K~ir + ir~ candidates is shown in Fig. 1 as points, and 
the open histogram is the distribution for the simulated 
K + K~tt + it~ events. The simulated distribution is nor- 
malized to the data in the region Xkk; + - ^ ^0 where 
the backgrounds and radiative corrections are insignif- 
icant. The experimental distribution has contributions 
from background processes, but the simulated distribu- 
tion is also broader than the expected 3C % 2 distribution. 
This is due to multiple soft-photon emission from the ini- 
tial state and radiation from the final-state charged par- 
ticles, which are not taken into account by the fit, but are 
present in both data and simulation. The shape of the x 2 
distribution at high values was studied in detail [11, 12] 
using specific ISR processes for which a very clean sample 
can be obtained without any limit on the x 2 value. 

The cross-hatched histogram in Fig. 1 represents the 
background from e + e~ — ► qq events, which is based on 
the JETSET simulation. It is dominated by events with 
a hard ir° producing a fake ISR photon, and the similar 
kinematics cause it to peak at low values of Xkk + -• 



We evaluate this background in a number of E c . m . ranges 
by combining the ISR photon candidate with another 
photon candidate in both data and simulated events, 
and comparing the 7r° signals in the resulting 77 invari- 
ant mass distributions. The simulation gives an -E c . in- 
dependence consistent with the data, so we normalize 
it by an overall factor. The hatched histogram repre- 
sents the sum of this background and that from ISR 
e + e _ — » 7T + 7r _ 7r + 7r~ events with one or two misidentified 
7r ± , which also contributes at low x 2 values. We estimate 
the contribution as a function of E c , m from a simulation 
using the known cross section [11]. 

All remaining background sources are either negligible 
or give a Xkktt+tt- distribution that is nearly uniform 
over the range shown in Fig. 1. We therefore define a 
signal region X 2 KKir+-n-<^i and estimate the sum of the 
remaining backgrounds from the difference between the 
number of data and simulated entries in a control region, 
30 < XxKn+jr- < 60. This difference is normalized to the 
corresponding difference in the signal region, as described 
in detail in Refs. [11, 12]. The signal region contains 
34635 data and 14077 simulated events, and the control 
region contains 4634 data and 723 simulated events. 




FIG. 2: The invariant mass distribution for K K~n 7r _ can- 
didates in the data (points): the cross-hatched, hatched and 
open histograms represent, cumulatively, the non-ISR back- 
ground, the contribution from ISR ir + iv~ir + iv~ events, and 
the ISR background from the control region of Fig. 1. 

Figure 2 shows the K + K~Tt + Tr~ invariant mass dis- 
tribution from threshold up to 5.0 GeV/c 2 for events in 
the signal region. Narrow peaks are apparent at the J/ip 
and i/j(2S) masses. The cross-hatched histogram repre- 
sents the qq background, which is negligible at low mass 
but becomes large at higher masses. The hatched region 
represents the ISR 7r + 7r~7r + 7r~ contribution, which we es- 
timate to be 2.4% of the selected events on average. The 
open histogram represents the sum of all backgrounds, 
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including those estimated from the control region. They 
total 6-8% at low mass but account for 20-25% of the ob- 
served data near 4 GeV/c 2 and become the largest con- 
tribution near 5 GeV/c 2 . 

We subtract the sum of backgrounds in each mass bin 
to obtain a number of signal events. Considering un- 
certainties in the cross sections for the background pro- 
cesses, the normalization of events in the control region 
and the simulation statistics, we estimate a systematic 
uncertainty on the signal yield that is less than 3% in 
the 1.6-3 GeV/c 2 mass region, but increases to 3-5% in 
the region above 3 GeV/c 2 . 



B. Selection Efficiency 

The selection procedures applied to the data are also 
applied to the simulated signal samples. The resulting 
K + K~ir + ir~ invariant-mass distributions in the signal 
and control regions are shown in Fig. 3(a) for the phase 
space simulation. The broad, smooth mass distribution 
is chosen to facilitate the estimation of the efficiency as a 
function of mass, and this model reproduces the observed 
distributions of kaon and pion momenta and polar angles. 
We divide the number of reconstructed simulated events 
in each mass interval by the number generated in that 
interval to obtain the efficiency shown as the points in 
Fig. 3(b). The 3 rd order polynomial fit to the points is 
used for further calculations. We simulate events with 
the ISR photon confined to the angular range 20-160° 
with respect to the electron beam in the e + e~ cm. frame, 
which is about 30% wider than the EMC acceptance. 
This efficiency is for this fiducial region, but includes the 
acceptance for the final-state hadrons, the inefficiencies of 
the detector subsystems, and event loss due to additional 
soft-photon emission. 

The simulations including the 0(lO2O)7r + 7r~ and/or 
K + K~ fo(980) channels have very different mass and an- 
gular distributions in the K + K~ ir + ir~ rest frame. How- 
ever, the angular acceptance is quite uniform for ISR 
events, and the efficiencies are consistent with those from 
the phase space simulation within 3%. To study possible 
mis-modeling of the acceptance, we repeat the analysis 
with the tighter requirements that all charged tracks be 
within the DIRC acceptance, 0.45 < # c h < 2.4 radians, and 
the ISR photon be well away from the edges of the EMC, 
0.35 < 6>isr < 2.4 radians. The fraction of selected data 
events satisfying the tighter requirements differs from the 
simulated ratio by 3.7%. We conservatively take the sum 
in quadrature of this variation and the 3% model varia- 
tion (5% total) as a systematic uncertainty due to accep- 
tance and model dependence. 

We correct for mis-modeling of the shape of the 
Xkkh+tt- distribution by (3.0 ±2.0)% and the track find- 
ing efficiency following the procedures described in de- 
tail in Ref. [11]. We use a comparison of data and sim- 
ulated xlvr distributions in the much larger samples of 
ISR 7r + 7r~7r + 7r - events. We consider data and simulated 
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FIG. 3: (a) The invariant mass distributions for simulated 
K + K~ tt + tt~ events in the phase space model, reconstructed 
in the signal (open) and control (hatched) regions of Fig. 1; 
(b) net reconstruction and selection efficiency as a function 
of mass obtained from this simulation (the curve represents a 
3 rd order polynomial fit). 
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FIG. 4: The e + e~—>K + K~n + n~ cross section as a function 
of the effective e + e~ cm. energy measured with ISR data at 
BABAR (dots). The direct measurements from DM1 [13] are 
shown as the open circles. Only statistical errors are shown. 



events that contain a high-energy photon plus exactly 
three charged tracks and satisfy a set of kinematical cri- 
teria, including a good \ 2 from a kinematic fit under the 
hypothesis that there is exactly one missing track in the 
event. We find that the simulated track-finding efficiency 
is overestimated by (0.8 ± 0.5)% per track, so we apply 
a correction of +(3 ± 2)% to the signal yield. 
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TABLE I: Measurements of the e + e — » K + K tt + 7v cross section (errors are statistical only). 



£c. m . (GeV) a (nb) £ c . m . (GcV) a (nb) E c . m . (GeV) a (nb) E c . m . (GeV) a (nb) E c . m . (GeV) a (nb) 



1.4125 





00± 0. 


.02 


2. 


.1375 


2.83± 


.13 


2.8625 


0. 


.50+ 0.05 


3.5875 


0. 


.12+ 0.03 


4.3125 


0.04± 0.02 


1.4375 





.01±0. 


.02 


2. 


.1625 


2.71±0 


.12 


2.8875 


0. 


.51±0.05 


3.6125 


0. 


.13+0.03 


4.3375 


0.04± 0.02 


1.4625 





.00± 0. 


.02 


2 


.1875 


2.46± 


.12 


2.9125 


0. 


.54± 0.05 


3.6375 


0. 


.12+0.03 


4.3625 


0.03± 0.02 


1.4875 





.04± 0. 


.02 


2 


.2125 


1.84± 


.10 


2.9375 


0. 


.46+ 0.05 


3.6625 


0. 


.11+0.03 


4.3875 


0.06± 0.02 


1.5125 





.03± 0. 


02 


2 


.2375 


1.66± 


.10 


2.9625 


0. 


45+ 0.05 


3.6875 


0. 


.25+ 0.03 


4.4125 


0.01± 0.02 


1.5375 





.11± 0. 


03 


2 


.2625 


1.59+0 


.09 


2.9875 


0. 


.46+ 0.05 


3.7125 


0. 


07+ 0.03 


4.4375 


0.03± 0.02 


1.5625 





,15± 0. 


.04 


2 


.2875 


1.66± 


.09 


3.0125 


0. 


.36+ 0.04 


3.7375 


0. 


.08+ 0.02 


4.4625 


0.06± 0.02 


1.5875 





.32±0. 


.05 


2 


.3125 


1.50+0 


.09 


3.0375 


0. 


.39+ 0.04 


3.7625 


0. 


.11+0.03 


4.4875 


0.03±0.02 


1.6125 





,48± 0. 


.06 


2 


.3375 


1.65+0 


.09 


3.0625 


0. 


.31+0.04 


3.7875 


0. 


.11+0.03 


4.5125 


0.04± 0.02 


1.6375 





.85± 0. 


.08 


2 


.3625 


1.56± 


.09 


3.0875 


2. 


.95+ 0.10 


3.8125 


0. 


.10+ 0.03 


4.5375 


0.01± 0.02 


1.6625 


1. 


,42± 0. 


.10 


2 


.3875 


1.49± 


.09 


3.1125 


1. 


.51+0.08 


3.8375 


0. 


.08+ 0.02 


4.5625 


0.02+ 0.02 


1.6875 


1 


.86± 0. 


.11 


2 


.4125 


1.46± 


.09 


3.1375 


0. 


.37+ 0.04 


3.8625 


0. 


.12+ 0.03 


4.5875 


0.05± 0.02 


1.7125 


2 


.36+0. 


.13 


2 


.4375 


1.48± 


.09 


3.1625 


0. 


.35+ 0.04 


3.8875 


0. 


.09+ 0.02 


4.6125 


0.02± 0.02 


1.7375 


2 


,67± 0. 


.13 


2 


.4625 


1.17±0 


.08 


3.1875 


0. 


.28± 0.04 


3.9125 


0. 


.09+ 0.02 


4.6375 


0.01±0.02 


1.7625 


3 


.51±0. 


15 


2 


.4875 


1.16+0 


.08 


3.2125 


0. 


.35+ 0.04 


3.9375 


0. 


.08+ 0.02 


4.6625 


0.04+ 0.02 


1.7875 


3 


,98± 0. 


.16 


2 


.5125 


1.21± 


.08 


3.2375 


0. 


31+0.04 


3.9625 


0. 


.10+0.02 


4.6875 


0.02± 0.02 


1.8125 


4 


,10± 0. 


.16 


2 


.5375 


0.94± 


.07 


3.2625 


0. 


.30+ 0.04 


3.9875 


0. 


.04+ 0.02 


4.7125 


0.03+ 0.02 


1.8375 


4. 


,68± 0. 


.17 


2 


.5625 


0.95±0 


.07 


3.2875 


0. 


.24+ 0.04 


4.0125 


0. 


.06+ 0.02 


4.7375 


0.01±0.02 


1.8625 


4. 


,49± 0. 


.17 


2 


.5875 


0.84± 


.07 


3.3125 


0. 


.22+ 0.04 


4.0375 


0. 


.07+ 0.02 


4.7625 


0.02± 0.02 


1.8875 


4 


.26±0. 


.17 


2 


.6125 


0.85±0 


.07 


3.3375 


0. 


.25+ 0.04 


4.0625 


0. 


.05+0.02 


4.7875 


0.01±0.02 


1.9125 


4. 


.30± 0. 


.16 


2 


.6375 


0.90± 


.07 


3.3625 


0. 


.16+0.03 


4.0875 


0. 


.06+ 0.02 


4.8125 


0.00± 0.02 


1.9375 


4. 


,20±0. 


.16 


2 


.6625 


0.82±0 


.06 


3.3875 


0. 


.17+0.03 


4.1125 


0. 


.06+ 0.02 


4.8375 


0.02± 0.02 


1.9625 


4. 


.13+0. 


.16 


2 


.6875 


0.70± 


.06 


3.4125 


0. 


.18+0.03 


4.1375 


0. 


.05+0.02 


4.8625 


0.00+ 0.02 


1.9875 


3 


,74± 0. 


.15 


2 


.7125 


0.86± 


.06 


3.4375 


0. 


.12+0.03 


4.1625 


0. 


06+ 0.02 


4.8875 


0.04± 0.02 


2.0125 


3 


,45± 0. 


.15 


2 


.7375 


0.81±0 


.06 


3.4625 


0. 


.17+0.03 


4.1875 


0. 


.05+0.02 


4.9125 


0.05±0.02 


2.0375 


3 


,38± 0. 


.14 


2 


.7625 


0.76± 


.06 


3.4875 


0. 


.17+0.03 


4.2125 


0. 


05+0.02 


4.9375 


0.02±0.02 


2.0625 


3 


,17±0. 


.14 


2 


.7875 


0.73± 


.06 


3.5125 


0. 


21+0.03 


4.2375 


0. 


.08+ 0.02 


4.9625 


0.00± 0.02 


2.0875 


3 


.23±0. 


14 


2 


.8125 


0.64± 


.05 


3.5375 


0. 


.14+0.03 


4.2625 


0. 


.04+ 0.02 


4.9875 


0.04± 0.02 


2.1125 


3 


,15±0. 


.14 


2 


.8375 


0.56±0 


.05 


3.5625 


0. 


.16+0.03 


4.2875 


0. 


.08+ 0.02 







We correct the simulated kaon identification efficiency 
using e + e~^> </>(1020)7^ K + K~-y events. Events with a 
hard ISR photon and two charged tracks, one of which is 
identified as a kaon, with a K + K~ invariant mass near 
the 4> mass provide a very clean sample, and we com- 
pare the fractions of data and simulated events with the 
other track also identified as a kaon, as a function of mo- 
mentum. The data-simulation efficiency ratio averages 
0.990 ± 0.001 in the 1-5 GeV/c momentum range with 
variations at the 0.01 level. We conservatively apply a 
correction of +(1.0 ± 1.0)% per kaon, or +(2.0 ± 2.0)% 
to the signal yield. 



C. Cross Section for e + e 



K+K-TT + TT- 



We calculate the e+e -^K + K tt + tt cross section as 
a function of the effective cm. energy from 



VKKir+ir- (E c .m. ) — 



dN 



KKn+-rt-~f i.E c .va 



dC(E c . m .) (E c .m.) 



(2) 



where i? c . m . = mmj+jt-c 2 , rnKKir+n- is the measured 
invariant mass of the K + K~tt + tt~ system, dN KKv + v - 7 
is the number of selected events after background sub- 
traction in the interval dE c , m ,, and exxi+jr-^c.m.) is 



the corrected detection efficiency. We calculate the differ- 
ential luminosity, d£(E c , m ,), in each interval dE c , m , from 
ISR /i + /x~7 events with the photon in the same fiducial 
range used for the simulation; the procedure is described 
in Refs. [11, 12]. From data-simulation comparison we 
conservatively estimate a systematic uncertainty on dC 
of 3%. This dC has been corrected for vacuum polar- 



TABLE II: Summary of corrections and systematic uncer- 
tainties on the e + e~^ K + K~-k + -k~ cross section. The total 
correction is the linear sum of the components and the total 
uncertainty is the sum in quadrature. 



Source Correction Uncertainty 


Rad. Corrections 




1% 






Backgrounds 




3%, 


m K Ktt+tt- 


< 3 GeV/c 2 






5%, 


m KKlT+TT- 


> 3 GeV/c 2 


Model Dependence 




5% 






Xkk,+ t - Distn - 


+3% 


2% 






Tracking Efficiency 


+3% 


2% 






Kaon ID Efficiency 


+2% 


2% 






ISR Luminosity 




3% 






Total 


+8% 


7%, 
9%, 


m KKlT+TT- 

m KKir+ir- 


< 3 GeV/c 2 
> 3 GeV/c 2 
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ization and final-state soft-photon emission; the former 
should be excluded when using these data in calculations 
of 57,-2. 

For the cross section measurement we use the tighter 
angular criteria on the charged tracks and the ISR pho- 
ton, discussed in Sec. IV B, to exclude possible errors 
from incorrect simulation of the EMC and DCH edge ef- 
fects. We show the cross section as a function of E c . m . in 
Fig. 4, with statistical errors only, and provide a list of 
our results in Table I. The result is consistent with the 
direct measurement by DM1 [13], and with our previous 
measurement of this channel [11] but has much better 
statistical precision. The systematic uncertainties, sum- 
marized in Table II, affect the normalization, but have 
little effect on the energy dependence. 

The cross section rises from threshold to a peak value 
of about 4.7 nb near 1.85 GeV, then generally decreases 
with increasing energy. In addition to narrow peaks at 
the J/ip and ip(2S) masses, there are several possible 
wider structures in the 1.8-2.8 GeV region. Such struc- 
tures might be due to thresholds for intermediate reso- 
nant states, such as </>/o(980) near 2 GeV. Gaussian fits 
to the simulated line shapes give a resolution on the mea- 
sured K + K~ir + ir~ mass that varies between 4.2 MeV/c 2 
in the 1.5-2.5 GeV/c 2 region and 5.5 MeV/c 2 in the 2.5- 
3.5 GeV/c 2 region. The resolution function is not purely 
Gaussian due to soft-photon radiation, but less than 10% 
of the signal is outside the 25 MeV/c 2 mass bin. Since 
the cross section has no sharp structure other than the 
J/ij) and if}{2S) peaks discussed in Sec. VIII below, we 
apply no correction for resolution. 
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FIG. 5: (a) Invariant mass of the K~n + pair versus that 
of the K + ty~ pair; (b) sum of projections of (a); (c) sum of 
projections of the if* bands of (a), with events in the overlap 
region taken only once. The (f>ir + Tv~ and KKp are vetoed. 



in Fig. 6(b). The bulk of Fig. 6(b) shows a strong posi- 
tive correlation, characteristic of K*°Kir final states with 
no higher resonances. The horizontal band in Fig. 6(b) 
corresponds to the peak region in Fig. 6(a), and is consis- 
tent with contributions from the if i (1270) and If i (1400) 
resonances. There is also an indication of a vertical band 
in Fig. 6(b), perhaps corresponding to a K*°K resonance 
at -1.5 GeV/c 2 . 



D. Substructure in the K + K~ir + n~ Final State 



Our previous study [11] showed many intermediate res- 
onances in the K + K ~tt + it~ final state. With the larger 
data sample used here, they can be seen more clearly 
and, in some cases, studied in detail. Figure 5(a) shows 
a scatter plot of the invariant mass of the K~n + pair 
versus that of the K + ir~ pair, and Fig. 5(b) shows the 
sum of the two projections. Here we have suppressed the 
contributions from 4>n + Tr~ and if + K ~p(770) by requir- 
ing \m(K+K-)-m(4>)\ > 10 MeV/c 2 and ^(tt+tt") - 
m(p)\ > 100 MeV/c 2 , where m{4>) and m(p) values are 
taken from the Particle Data Group (PDG) tables [5]. 
Bands and peaks corresponding to the if* (892) and 
if*°(1430) are visible. In Fig. 5(c) we show the sum 
of projections of the if *°(892) bands, defined by lines in 
Fig. 5(a), with events in the overlap region plotted only 
once. No if* (892) signal is seen, confirming that the 
e+e _ ->if* (892)F*°(892) cross section is small. We ob- 
serve associated if*°(892)if *,°(1430) production, but it 
is mostly from J/ip decays (see Sec. VIII). 

We combine K /K*° candidates within the lines in 
Fig. 5(a) with the remaining pion and kaon to ob- 
tain the if *°7r H invariant mass distribution shown in 
Fig. 6(a), and the if * 7r+" vs. K*°K~ + mass scatter plot 




1 1.5 2 2.5 . 
m(KV") (GeV/c 2 ) 



1 1.5 2 2.5 3 

m(K*K + ") (GeV/c 2 ) 



FIG. 6: (a) The if *°tt invariant mass distribution; (b) the 
if *°7r mass versus K*°K mass. 

We now suppress K*°Kir by considering only events 
outside the lines in Fig. 5(a). In Fig. 7 we show the 
K^'k+tt^ invariant mass (two entries per event) vs. that 
of the 7r + 7r~ pair, along with its two projections. There 
is a strong p(770) -^7r + 7r~ signal, and the if ± 7r + 7r~ mass 
projection shows further indications of the if i (1270) and 
ifi(1400) resonances, both of which decay into ifp(770). 
There are suggestions of additional structure in the 7r + 7r~ 
mass distribution, including possible /o(980) shoulder 
and a possible enhancement near the /2(1270), however 
the current statistics do not allow us to make definitive 
statements. 

The separation of all these, and any other, intermedi- 
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m(7t + 7t ) (GeV/c ) 

FIG. 7: (a) Invariant mass of the K ± it + tt~ combinations ver- 
sus that of the 7r + 7r~ pair; (b) the K ± tt + -k~ and (c) 7r + 7r~ 
mass projections of (a). 

ate states involving relatively wide resonances requires a 
partial wave analysis. This is beyond the scope of this 
paper. Here we present the cross section for the sum of 
all states including a if* (892), and study intermediate 
states that include a narrow <p or /q resonance. 

E. The e + e-—yK*°KTT Cross Section 

Signals for the if* (892) and A 2 *°(1430) are clearly vis- 
ible in the K if* mass distributions in Fig. 5(b) and, 
with a different bin size, in Fig. 8(a). We perform a 
fit to this distribution using P-wave Breit-Wigner (BW) 
functions for the K*° and if| signals and a third-order 
polynomial function for the remainder of the distribu- 
tion taking into account the Kir threshold. The result 
is shown in Fig. 8(a). The fit yields a K*° signal of 
19738 ± 266 events with m(Kn) = 896.2 ± 0.3 MeV/c 2 
and r(ifvr) = 50.6 ± 0.9 MeV, and a A 2 *° signal of 
1786± 127 events with m(Kn) = 1428.5±3.9 MeV/c 2 and 
r(Kn) = 113.7 ± 9.2 MeV. These values are consistent 
with current world averages [5], and the fit describes the 
data well, indicating that contributions from any other 
resonances decaying into K ± 7r =F are small. 

We perform a similar fit to the data in bins of the 
K + K ~7r + 7r~ invariant mass, with the resonance masses 
and widths fixed to the values obtained by the overall fit. 
Since there is at most one K*° per event, we convert the 
resulting K*° yield in each bin into an "inclusive" e + e~^ 
K*°Kir cross section, following the procedure described 
in Sec. IV C. This cross section is shown in Fig. 8(b) and 
listed in Table III for the effective cm. energies from 
threshold up to 3.5 GeV. At higher energies the signals 
are small and contain an unknown, but possibly large, 
contribution from e + e~^qq events. There is a rapid rise 
from threshold to a peak value of about 4 nb at 1.84 GeV, 




E c . m . (GeV) 

FIG. 8: (a) The K ± ty^ 1 mass distribution (two entries per 
event) for all selected K + K~ -k + tt~ events: the solid line rep- 
resents a fit including two resonances and a polynomial func- 
tion (see text), shown separately as the dashed line; (b) the 
cross section obtained from the A* (892) sig- 
nal by a similar fit in each 25 MeV/c 2 mass bin. 

followed by a very rapid decrease with increasing energy. 
There are suggestions of narrow structure in the peak 
region, but the only statistically significant structure is 
the J I if) peak, which is discussed below. 

The e + e~^ K* a Kir contribution is a large fraction of 
the total K + K~tt + it~ cross section at all energies above 
its threshold, and dominates in the 1.8-2.0 GeV region. 
We are unable to extract a meaningful measurement of 
the K%°Kir cross section from this data sample because 
it is more than ten times smaller. The K + K~ p°(770) in- 
termediate state makes up the majority of the remainder 
of the cross section and it can be estimated as a difference 
of the values in Tabic I and Table III for the K + K~ir + ir~ 
and K*°Kir final states. 



F. The 0(lO2O)7r+7r~ Intermediate State 

Intermediate states containing relatively narrow reso- 
nances can be studied more easily. Figure 9(a) shows 
a scatter plot of the invariant mass of the ir + ir~ pair 
versus that of the K + K~ pair. Horizontal and vertical 
bands corresponding to the p°(770) and <f>, respectively, 
are visible, and there is a concentration of entries on the 
<f> band corresponding to the correlated production of <f> 
and /o(980). The <j> signal is also visible in the K + K~ 
mass projection, Fig. 9(c). The large contribution from 
the /?(770), coming from the K\ decay, is nearly uniform 
in the K + K~ mass, and the cross-hatched histogram 
shows the non-A + A~7r + 7r~ background estimated from 
the control region in XfcKn+n- • The coss-hatched his- 
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TABLE III: Measurements of the e+e" -> K*°{892)K n cross section (errors are statistical only). 
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togram also shows a <j) peak, but this is a small fraction 
of the events. Subtracting this background and fitting 
the remaining data gives 1706±56 events produced via 
the </>7r + 7r~ intermediate state. 

To study the 07r + 7r~ channel, we select candidate 
events with a K + K~ invariant mass within 10 MeV/c 2 
of the <f> mass, indicated by the inner vertical lines in 
Figs. 9(a,c), and estimate the non-0 contribution from 
the mass sidebands between the inner and outer verti- 
cal lines. In Fig. 9(b) we show the 7r + 7r~ invariant mass 
distributions for <f> candidate events, sideband events and 
X 2 control region events as the open, hatched and cross- 
hatched histograms, respectively, and in Fig. 9(d) we 
show the numbers of entries from the candidate events 
minus those from the sideband and control region. There 
is a clear /q peak over a broad mass distribution, with 
no indication of associated p° production. 

A coherent sum of two Breit-Wigner functions is suf- 
ficient to describe the invariant mass distribution of the 
7r + 7r~ pair recoiling against a <\> in Fig. 9(d). We fit with 
the function: 

F(m) 
Ai(m) 



y/1 -4m 2 /to 2 • |Ai(to) +e^A 2 (m)| 2 , (3) 
rriiTi \f~N~i/ (to 2 - to 2 + irriiTi) , 



where m is the 7r + 7r~ invariant mass, to^ and Ti are the 
parameters of the i th resonance, tp is their relative phase 
and Ni are normalization parameters, corresponding to 
the number of events under each BW. One BW corre- 
sponds to the /o(980), but a wide range of values of the 
other parameters can describe the data. Fixing the rela- 
tive phase to ip — tt and the parameters of the first BW 
to toi = 0.6 GeV/c 2 and T x = 0.45 GeV (which could be 
interpreted as the /o(600) [5]), we obtain the fit shown in 
Fig. 9(d). It describes the data well and gives an /o(980) 
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FIG. 9: (a) The it tt~ vs. the K K~ invariant masses 
for all selected K + K~ -k + tt~ events; (b) the -k + -k~ invari- 
ant mass projections for events in the <f> peak (open his- 
togram), sidebands (hatched) and background control region 
(cross-hatched); (c) the K + K~ mass projections for all events 
(open) and control region (cross-hatched); (d) the difference 
between the open and the sum of the other histograms in (b). 



signal of 262±30 events, with m 2 = 0.973±0.003 GeV/c 2 
and T 2 = 0.065 ± 0.013 GeV, consistent with the PDG 
values [5]. There is a suggestion of an / 2 (1270) peak in 
the data, but it is much smaller than the /o peak and we 
do not consider it further. 

We obtain the number of e + e~ — > (/>7r + 7r~ events in 
bins of <pTT + TT~ invariant mass by fitting the K + K~ in- 
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variant mass projection in that bin after subtracting non- 
K + K~tt + it~ background. Each projection is a subset of 
Fig. 9(c), where the curve represent a fit to the full sam- 
ple. In each mass bin, all parameters are fixed to the 
values obtained from the overall fit except the numbers 
of events in the <j> peak and the non-</> component. 

The efficiency may depend on the details of the pro- 
duction mechanism. Using the two-pion mass distribu- 
tion in Fig. 9(d) as input, we simulate the tt + it~ system 
as an S-wave comprising two scalar resonances, with pa- 
rameters set to the values given above. To describe the 
4>it + it~ mass distribution we use a simple model with one 
resonance, the </>(1680), of mass 1.68 GeV/c 2 and width 
0.2 GeV, decaying to (j)f . The simulated reconstructed 
spectrum is shown in Fig. 10(a). There is a sharp in- 
crease at about 2 GeV/c 2 due to the 0/o(98O) threshold. 
All other structure is determined by phase space and a 
to~ 2 falloff with increasing mass. 



,4000 




m(KTCW) (GeV/c ) 



FIG. 10: (a) The K K~n n~ invariant mass distributions 
from the 07r + 7r~ simulation described in the text, recon- 
structed in the signal (open) and control (hatched) regions; 
(b) net reconstruction and selection efficiency as a function 
of mass: the solid line represents a cubic fit, and the dashed 
line the corresponding fit for the space phase model shown in 
Fig. 3. 



Dividing the number of reconstructed events in each 
bin by the number of generated ones, we obtain the effi- 
ciency as a function of 07r + 7r~ mass shown in Fig. 10(b). 
The solid line represents a fit to a third order polyno- 
mial, and the dashed line the corresponding fit to the 
phase space model from Fig. 3. The model dependence 
is weak, giving confidence in the efficiency calculation. 
We calculate the e + e~ — > 07r + 7r~ cross section as de- 
scribed in Sec. IV C but using the efficiency from the fit 
to Fig. 10(b) and dividing by the <f)—>K + K~ branching 
fraction of 0.491 [5]. We show our results as a func- 
tion of energy in Fig. 11 and list them in Table IV. 



The cross section has a peak value of about 0.6 nb at 
about 1.7 GeV, then decreases with increasing energy 
until 0(1O2O)/ O (98O) threshold, around 2.0 GeV. From 
this point it rises, falls sharply at about 2.2 GeV, and 
then decreases slowly. Except in the charmonium region, 
the results at energies above 2.9 GeV are not meaningful 
due to small signals and potentially large backgrounds, 
and are omitted from Table IV. Figure 11 displays the 
cross-section up to 4.5 GeV to show the signals from the 
J/ip and ^>(2S) decays. They are discussed in Sec. VIII. 
There are no previous measurements of this cross section, 
and our results are consistent with the upper limits given 
in Ref. [13]. 




E c . m . (GeV) 



FIG. 11: The e + e — > <fnx + -K cross section as a function of 
the effective e + e~ cm. energy. 



We perform a study of the angular distributions in the 
</i>(1020)7r + 7r~ final state by considering all K + K~ii + 'k~ 
candidate events with mass below 3 GeV/c 2 , binning 
them in terms of the cosine of the angles defined below, 
and fitting the background-subtracted K + K~ mass pro- 
jections. The efficiency is nearly uniform in these angles, 
so we study the number of events in each bin. We define 
the <j) production angle, as the angle between the <j> 
momentum and the e~ beam direction in the rest frame 
of the (f>n + ir~ system. The distribution of cosO^, shown 
in Fig. 12(a), is consistent with the uniform distribution 
expected if S-wave two-body channels <pX, X — > 7r + 7r~ 
dominate the <f>ir + iT~ system. We define the pion and 
kaon helicity angles, 6 W + and &k+ as those between the 
7T+ and the 7r + 7r~-system momenta in the tt + tt~ rest 
frame and between the K + and ISR photon momenta 
in the <j> res t frame, respectively. The distributions of 
cos 9^+ and cos 9^+, shown in Figs. 12(b) and 12(c), re- 
spectively, are consistent with those expected from scalar 
and vector meson decays. 
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TABLE IV: Measurements of the e + e — > <7j(1020)7r + 7r cross section (errors are statistical only). 
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FIG. 12: Distributions of the cosines of (a) the <f> production angle, (b) the pion helicity angle, and (c) the kaon helicity angle 
(see text) for e + e~ — > (f>ir + tv~ events: the lines represent the distributions expected if the it + tt~ system recoiling against a vector 
4> meson is produced in an S-wave, normalized to the number of events in the data. 



G. The 0(1O2O)/ O (98O) Intermediate State 

The narrow / (980) peak seen in Fig. 9(d) allows the 
selection of a fairly clean sample of (j)f events. We re- 
peat the analysis just described with the additional re- 
quirement that the 7r + 7r~ invariant mass be in the range 
0.85-1.10 GeV/c 2 . The fit to the full sample yields about 
700 events; all of these contain a true 0, but about 10% 
are from e+e~ — > </>7r + 7r~ events where the pion pair is 
not produced through the /o(980). 

We convert the numbers of fitted events in each mass 
bin into a measurement of the e+e~ — > </>(1020).fo(980) 
cross section as described above and dividing by the 
/o— > n + n~ branching fraction of two-thirds. The cross 
section is shown in Fig. 13 as a function of the effective 
cm. energy and is listed in Table V. Its behavior near 
threshold docs not appear to be smooth, but is more 
consistent with a steep rise to a value of about 0.3 nb at 
1.95 GeV followed by a slow decrease that is interrupted 
by a structure around 2.175 GeV. Possible interpreta- 
tions of this structure are discussed in Sec. VII. Again, 



the values are not meaningful for the effective cm. above 
about 2.9 GcV/c 2 , except for the J/tp and ip(2S) signals, 
discussed in Sec. VIII. 



V. THE K+K~ ttV FINAL STATE 
A. Final Selection and Backgrounds 

The K + K~ir°ir° sample contains background from 
the ISR processes e + e~ — > K + K~Tr°-f and K + K~rjj, 
in which two soft photon candidates from machine- or 
detector-related background combine with the relatively 
energetic photons from the 7r° or 77 to form two fake ir° 
candidates. We reduce this background using the helic- 
ity angle between each reconstructed ir° direction and the 
direction of its higher-energy photon daughter calculated 
in its rest frame. If the cosines of both helicity angles are 
higher than 0.85, we remove the event. 

Figure 14 shows the distribution of Xkkk * ^ or 
the remaining candidates together with the simulated 
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TABLE V: Summary of the e + e — > 0(lO2O)/o(98O) cross section measurement. Errors are statistical only. 
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a (nb) 


1.8875 


0.02 ± 0.02 


2.1625 


0.56 ± 0.08 


2.4375 


0.13 ± 0.04 


2.7125 


0.05 ± 0.02 


1.9125 


0.07 ± 0.04 


2.1875 


0.38 ± 0.07 


2.4625 


0.09 ± 0.04 


2.7375 


0.02 ± 0.01 


1.9375 


0.30 ± 0.07 


2.2125 


0.21 ± 0.06 


2.4875 


0.08 ± 0.03 


2.7625 


0.06 ± 0.02 


1.9625 


0.27 ± 0.07 


2.2375 


0.17 ± 0.05 


2.5125 


0.09 ± 0.03 


2.7875 


0.03 ± 0.02 


1.9875 


0.25 ± 0.07 


2.2625 


0.12 ± 0.04 


2.5375 


0.07 ± 0.03 


2.8125 


0.04 ± 0.02 


2.0125 


0.25 ± 0.07 


2.2875 


0.13 ± 0.05 


2.5625 


0.02 ± 0.02 


2.8375 


0.05 ± 0.02 


2.0375 


0.28 ± 0.07 


2.3125 


0.12 ± 0.04 


2.5875 


0.06 ± 0.03 


2.8625 


0.02 ± 0.02 


2.0625 


0.43 ± 0.08 


2.3375 


0.14 ± 0.04 


2.6125 


0.07 ± 0.03 


2.8875 


0.01 ± 0.01 


2.0875 


0.28 ± 0.07 


2.3625 


0.13 ± 0.05 


2.6375 


0.06 ± 0.03 


2.9125 


0.02 ± 0.01 


2.1125 


0.54 ± 0.09 


2.3875 


0.13 ± 0.04 


2.6625 


0.05 ± 0.03 


2.9375 


0.00 ± 0.00 


2.1375 


0.46 ± 0.08 


2.4125 


0.14 ± 0.05 


2.6875 


0.03 ± 0.02 


2.9625 


0.02 ± 0.01 
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1.5 2 2.5 3 3.5 4 

E c . m . (GeV) 

FIG. 13: The e+e"^ 0(1O2O)/ O (98O) cross section as a func- 
tion of the effective e + e~ cm. energy obtained from the 
K + K~ tt + -k~ final state. 

K + K~ir°ir° events. Again, the distributions are broader 
than those for a typical 6C % 2 due to higher order ISR, 
and we normalize the histogram to the data in the region 
Xkk^tt < 10- The cross-hatched histogram in Fig. 14 
represents background from e + e~-^qq events, evaluated 
in the same way as for the K + K~tt + tt~ final state. The 
hatched histogram represents the sum of this background 
and that from ISR 7r + 7r~7r°7r° events with both charged 
pions misidentified as kaons, evaluated using the simula- 
tion. 

The dominant background in this case is from resid- 
ual ISR K+R-tt and K + K~-q events, as well as ISR 
K + K~ir°ir°ir events. Their simulated contribution, 
shown as the dashed histogram in Fig. 14, is consistent 
with the data in the high Xkkw°tt° region. All other back- 
grounds are either negligible or distributed uniformly in 

Xkkk°-k°- We denne a si S nal region, Xkk-k°w < 50, con- 
taining 4425 data and 6948 simulated events, and a con- 
trol region, 50 < X'kkw -* < containing 1751 data 




FIG. 14: Distribution of x 2 from the six-constraint fit for 
K + K~ tt°-k° candidates in the data (points). The open his- 
togram is the distribution for simulated signal events, nor- 
malized as described in the text. The cross-hatched, hatched 
and dashed histograms represent, cumulatively, the back- 
grounds from non-ISR events, ISR tt + tt~ -k°tv° events, and ISR 
K+K-tv , K+K-j) and K+K "ttW events. 



and 848 simulated events. 

Figure 15 shows the K + K~it it invariant mass dis- 
tribution from threshold up to 5 GeV/c 2 for events in 
the signal region. The qq background (cross-hatched his- 
togram) is negligible at low masses but forms a large frac- 
tion of the selected events above about 4 GeV/c 2 . The 
ISR 7r + 7r _ 7r°7r° contribution (hatched region) is negligi- 
ble except in the 1.5-2.5 GeV/c 2 region. The sum of all 
other backgrounds, estimated from the control region, 
is the dominant contribution below 1.6 GeV/c 2 and non 
negligible everywhere. The total background in the 1.6- 
2.5 GeV/c 2 region is 15-20% (open histogram). 

We subtract the sum of backgrounds from the number 
of selected events in each mass bin to obtain a number 
of signal events. Considering uncertainties in the cross 
sections for the background processes, the normalization 
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m(K + K~JlV) (GeV/c 2 ) 



7r°7r° candi- 



FIG. 15: Invariant mass distribution for K K~ 
dates in the data (points). The cross-hatched, hatched and 
open histograms represent, cumulatively, the non-ISR back- 
ground, the contribution from ISR 7r + 7r~7r°7r° events, and the 
ISR background from the control region of Fig. 14. 
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FIG. 16: (a) Invariant mass distribution for simulated 
K + if ~7r%° events in the signal (open) and control (hatched) 
regions (see Fig. 14); (b) net reconstruction and selection ef- 
ficiency as a function of mass obtained from this simulation 
(the curve represents a third order polynomial fit). 



of events in the control region and the simulation statis- 
tics, we estimate a systematic uncertainty on the signal 
yield after background subtraction as less than 5% in the 
1.6-3.0 GeV/c 2 region, but increases to 10% in the region 
above 3 GeV/c 2 . 



B. Selection Efficiency 

The detection efficiency is determined in the same 
manner as in Sec. IV B. Figure 16(a) shows the sim- 
ulated K + K~Tr°Tr a invariant mass distributions in the 
signal and control regions from the phase space model. 
We divide the number of reconstructed events in each 
40 MeV/c 2 mass interval by the number generated ones 
in that interval to obtain the efficiency shown as the 
points in Fig. 16(b); a third order polynomial fit to the 
efficiency is used to calculate the cross section. Again, 
the simulation of the ISR photon covers a limited angu- 
lar range, about 30% wider than EMC acceptance, and 
shown efficiency is factor 0.7 lower than for the hadronic 
system alone. Simulations assuming dominance of the 
4>^K + K~ and/or fo-^>ir°Tr a channels give consistent re- 
sults, and we apply the same 5% systematic uncertainty 
for possible model dependence as in Sec. IV B. 

We correct for mis-modeling of the track finding and 
kaon identification efficiencies as in Sec. IV B, and for 
the shape of the Xkktt tt° distribution analogously, using 
the result in Ref. [12], (0 ± 6)%. We correct the tt°- 
finding efficiency using the procedure described in detail 
in Ref . [12]. From ISR e + e~^ii>7r 7^7r + 7r~7r 7r°7 events 
selected with and without the ir° from the u> decay, we 



find that the simulated efficiency for one 7r° is too high 
by (2.8±1.4)%. Conservatively we apply a correction of 
+ (5.6 ± 2.8)% for two tt° in the event. 



C. Cross Section for e+e" -> K + K-it°Tr° 

We calculate the cross section for e+e~ — > K + K ~7r°7r° 
in 40 McV E c . m . intervals from the analog of Eq. 2, using 
the invariant mass of the K + K ~7r°7r° system to deter- 
mine the effective cm. energy. We show the first mea- 
surement of this cross section in Fig. 17 and list the re- 
sults obtained in Table VI. The cross section rises to a 
peak value near 1 nb at 2 GeV, falls sharply at 2.2 GeV, 
then decreases slowly. The only statistically significant 
structure is the J/tp peak. The drop at 2.2 GeV is sim- 
ilar to that seen in the K + K~tt + tt~ mode. Again, dC 
includes corrections for vacuum polarization that should 
be omitted from calculations of g M — 2. 

The simulated K + K ~7r°7r° invariant mass resolution 
is 8.8 MeV/c 2 in the 1.5-2.5 GeV/c 2 mass range, and in- 
creases with mass to 11.2 MeV/c 2 in the 2.5-3.5 GeV/c 2 
range. Since less than 20% of the events in a 40 MeV/c 2 
bin are reconstructed outside that bin and the cross sec- 
tion has no sharp structure other than the J/tp peak, we 
again make no correction for resolution. The point-to- 
point systematic errors are much smaller than statistical 
ones, and the errors on the normalization are summa- 
rized in Table VII, along with the corrections that were 
applied to the measurements. The total correction is 
+9.2%, and the total systematic uncertainty is 10% at 
low mass, increasing to 14% above 3 GeV/c 2 . 
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TABLE VI: Measurements of the e + e — > K + K 7r°7r° cross section (errors are statistical only). 



E c . m . (GeV) 


a 


(nb) 


£ c .m. (GeV) 


a 


(nb) 


£cm. (GeV) 


a 


(nb) 


Bc.m. (GeV) 


o 


(nb) 


1.4200 


0.00 


± 


0.05 


2.3400 


0.35 


± 


0.06 


3.2600 


0.13 


± 


0.03 


4.1800 


0.02 


± 0.01 


1.4600 


0.12 


± 


0.07 


2.3800 


0.29 


± 


0.06 


3.3000 


0.09 


± 


0.03 


4.2200 


0.03 


± 0.01 


1.5000 


0.00 


± 


0.07 


2.4200 


0.38 


± 


0.06 


3.3400 


0.09 


± 


0.03 


4.2600 


0.03 


± 0.01 


1.5400 


0.01 


± 


0.08 


2.4600 


0.38 


± 


0.06 


3.3800 


0.08 


± 


0.02 


4.3000 


0.03 


± 0.01 


1.5800 


0.03 


± 


0.09 


2.5000 


0.22 


± 


0.05 


3.4200 


0.11 


± 


0.03 


4.3400 


0.03 


± 0.01 


1.6200 


0.09 


± 


0.09 


2.5400 


0.25 


± 


0.05 


3.4600 


0.06 


± 


0.02 


4.3800 


0.01 


± 0.01 


1.6600 


0.31 


± 


0.11 


2.5800 


0.25 


± 


0.05 


3.5000 


0.04 


± 


0.02 


4.4200 


0.05 


± 0.01 


1.7000 


0.35 


± 


0.11 


2.6200 


0.25 


± 


0.05 


3.5400 


0.06 


± 


0.02 


4.4600 


0.03 


± 0.01 


1.7400 


0.49 


± 


0.11 


2.6600 


0.28 


± 


0.05 


3.5800 


0.07 


± 


0.02 


4.5000 


0.04 


± 0.01 


1.7800 


0.51 


± 


0.12 


2.7000 


0.16 


± 


0.04 


3.6200 


0.04 


± 


0.02 


4.5400 


0.00 


± 0.01 


1.8200 


0.84 


± 


0.12 


2.7400 


0.22 


± 


0.04 


3.6600 


0.06 


± 


0.02 


4.5800 


0.02 


± 0.01 


1.8600 


0.94 


± 


0.11 


2.7800 


0.21 


± 


0.04 


3.7000 


0.08 


± 


0.02 


4.6200 


0.02 


± 0.01 


1.9000 


0.95 


± 


0.12 


2.8200 


0.13 


± 


0.04 


3.7400 


0.09 


± 


0.02 


4.6600 


0.02 


± 0.01 


1.9400 


0.80 


± 


0.11 


2.8600 


0.21 


± 


0.04 


3.7800 


0.02 


± 


0.02 


4.7000 


0.04 


± 0.01 


1.9800 


0.87 


± 


0.11 


2.9000 


0.11 


± 


0.03 


3.8200 


0.05 


± 


0.01 


4.7400 


0.02 


± 0.01 


2.0200 


1.00 


± 


0.10 


2.9400 


0.12 


± 


0.04 


3.8600 


0.04 


± 


0.01 


4.7800 


0.01 


± 0.01 


2.0600 


0.96 


± 


0.10 


2.9800 


0.12 


± 


0.04 


3.9000 


0.03 


± 


0.02 


4.8200 


0.01 


± 0.01 


2.1000 


0.90 


± 


0.10 


3.0200 


0.21 


± 


0.04 


3.9400 


0.02 


± 


0.01 


4.8600 


0.01 


± 0.01 


2.1400 


0.82 


± 


0.10 


3.0600 


0.16 


± 


0.04 


3.9800 


0.03 


± 


0.01 


4.9000 


0.03 


± 0.01 


2.1800 


0.58 


± 


0.08 


3.1000 


0.92 


± 


0.07 


4.0200 


0.05 


± 


0.01 


4.9400 


0.04 


± 0.02 


2.2200 


0.56 


± 


0.08 


3.1400 


0.19 


± 


0.04 


4.0600 


0.04 


± 


0.01 


4.9800 


0.04 


± 0.02 


2.2600 


0.37 


± 


0.07 


3.1800 


0.12 


± 


0.03 


4.1000 


0.03 


± 


0.01 








2.3000 


0.43 


± 


0.07 


3.2200 


0.14 


± 


0.03 


4.1400 


0.03 


± 


0.01 









TABLE VII: Summary of corrections and systematic uncer- 
tainties on the e + e~^ K + K~-K a n Q cross section. The total 
correction is the linear sum of the components and the total 
uncertainty is the sum in quadrature. 



Source 



Correction 



Uncertainty 



Rad. Corrections 




1% 




Backgrounds 




5%, m KK7r a 7V o 


< 3 GeV/c 2 
> 3 GeV/c 2 






10%, m KKn07r o 


Model Dependence 




5% 




Xh»o i0 Distn. 


0% 


6% 




Tracking Efficiency 


+1.6% 


0.8% 




Kaon ID Efficiency 


+2% 


2% 




7T° Efficiency 


+5.6% 


2.8% 




ISR Luminosity 




3% 




Total 


+9.2% 


10%, m KK ^aj> 


< 3 GeV/c 2 
> 3 GeV/c 2 






14%, m KKw07r o 



D. Substructure in the K+K~ tt°tt Final State 

A scatter plot of the invariant mass of the K~tt° versus 
that of the K + tt° pair is shown in Fig. 18(a) with two en- 
tries per event selected in the \ 2 signal region. Horizontal 
and vertical bands corresponding to the K* + (892) and 
K*~(892), respectively, are visible. Figure 18(b) shows 
as points the sum of the two projections of Fig. 18(a); 
a large K* ± (892) signal is evident. Fitting this distri- 
bution with the function discussed in Sec. IV E gives 
a good x 2 an d the curve shown on Fig. 18(b). The 
iT* ± (1430):i ; r ± (892) ratio is consistent with that for 




FIG. 17: The e+e"- 
of the effective e + e" 
BABAR. The errors are statistical only. 



(GeV) 



+ K K n tt cross section as a function 
cm. energy measured with ISR data at 



neutral K* seen in the K + K ~ir + ir~ channel, and the 
number of K* ± (892) in the peak is consistent with one 
per selected event. The hatched histogram in Fig. 18(b) 



represents the K 



i in events with the other K ^ir° 



mass within the lines in Fig. 18(a), but with events in the 
overlap region used only once, and shows no K* ± (892) 
signal. These results indicate that the e + e~^ K* ± K* Z ^ 
cross section is small and that the K* ± (892)K t tt° chan- 



19 





• 


1 






o 

> 

O 1.5 

O 








i (a) ; 


1 > 


























K"7t°) (GeV/c 2 ) 



FIG. 18: 

of the K 



1 1.5 2 
m(K + 7t°) (GeV/c 2 ) m(K T 7t u 

(a) Invariant mass of the K~tv° pair versus that 
7T° pair in selected K + K~ tv°ty° events (two entries 
per event); (b) sum of projections of (a) (dots, four entries per 
event). The curve represents the result of the fit described in 
the text. The hatched histogram is the K ± ty° distribution 
for events in which the other K T 7r° combination is within the 
if** (892) bands indicated in (a), with events in the overlap 
region taken only once. 



nels dominate the overall cross section. 

We find no signals for resonances in the K + K~ir° or 
K^tt decay modes. Since the K* ± (892)K^t: chan- 
nels dominate and the statistics are low in any mass bin, 
we do not attempt to extract a separate K* ± (892)K^n° 
cross section. The total K + K~ 7r°7r° cross section is 
roughly a factor of four lower than the if*°(892).R' ± 7r =F 
cross section observed in the K + K~ir + ir~ final state. 
This is consistent with what one might expect from 
isospin and the charged vs. neutral K* branching frac- 
tions into charged kaons. 



E. The 0(lO2O)7r°7r° Intermediate State 

The selection of events containing a 0(1020) — » K + K~ 
decay follows that in Section IV F. Figure 19(a) shows a 
scatter plot of the invariant mass of the 7r°7r° pair versus 
that of the K + K~ pair. A vertical band corresponding 
to the (j) is visible, whose intensity decreases with increas- 
ing 7r°7T mass except for an enhancement in the /o (980) 
region. The <fi signal is also visible in the K + K~ invari- 
ant mass projection shown in Fig. 19(c). The relative 
non-0 background is smaller than in the K + K ~ir + ir~ 
mode, but there is a large background from ISR 07r°, 
4>r] and/or 07r O 7r°7r° events, as indicated by the control 
region histogram (hatched) in Fig. 19(c). The contribu- 
tions from non-ISR and ISR 7r + 7r~7r°7r° events are negli- 
gible. Selecting cf) candidate and side band events as for 
the K + K~ir + n~ mode (vertical lines in Figs. 19(a,c)), 
we obtain the 7r°7r° mass projections shown as the open 
and cross-hatched histograms, respectively, in Fig. 19(b). 
Control region events (hatched histogram) are concen- 
trated at low masses. A peak corresponding to the 
/o(980) is visible over a relatively low background. 

In Fig. 19(d) we show the numbers of entries from 
the candidate events minus those from the sideband and 
control regions. A sum of two Breit-Wigner functions 
is again sufficient to describe the data. Fitting Eq. 3 
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FIG. 19: (a) Scatter plot of the it tv invariant mass vs. the 
K + K~ invariant mass for all selected K + K' events; (b) 

the 7T°7r invariant mass projections for events in the (j> peak 
(open histogram), sidebands (cross-hatched) and control re- 
gion (hatched); (c) the K + K~ mass projection for events in 
the signal (open) and control (hatched) regions; (d) difference 
between the open and other histograms in (b). 



with the parameters of one BW fixed to the values given 
in Sec. IV F, corresponding to the / (600), we obtain a 
good fit, shown as the curve in Fig. 19(d). This fit yields 
a /o(980) signal of 54±9 events with a mass m = 0.970± 
0.007 GeV/c 2 and width T = 0.081 ± 0.021 GeV consis- 
tent with PDG values [5] . Due to low statistics and high 
backgrounds, we do not extract an e + e~ — * 0(lO2O)7r°7r° 
cross section. 



F. The 0(1O2O)/ O (98O) Intermediate State 

Since the background under the / (980) peak in 
Figs. 19(b,d) is relatively low we are able to extract the 
0(1O2O)/ O (98O) contribution. As in Sec. IV G, we require 
the dipion mass to be in the range 0.85-1.10 GeV/c 2 and 
fit the background-subtracted K + K~ mass projection in 
each bin of K + K~ir°ir° mass to obtain a number of <pfo 
events. Again, about 10% of these are (/>7r 7r events in 
which the 7r°7r° pair is not produced through the /o, but 
this does not affect the conclusions. 

We convert the number of fitted events in each mass 
bin into a measurement of the e + e~ — > 0(lO2O)/o(98O) 
cross section as described above and dividing by the 
/o(980) — * TT°ir branching fraction of one-third. The 
cross section is shown in Fig. 20 as a function of E c . m , 
and is listed in Table VIII. Due to smaller number 
of events, we have used larger bins at higher energies. 
The overall shape is consistent with that obtained in the 
K + K~7t + it~ mode (see Fig. 13), and there is a sharp 
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TABLE VIII: Measurements of the e + e — * 0(1020) /o (980) cross section (/o — > 7r°7r°, errors are statistical only). 

gem. (GeV) n : nl.) E c . m . (GeV) a (nb) g c , m . (GeV) n :„|,:~ 

1.88-1.92 0.078+^*3 2.12-2.16 0.397+^ 2.44-2.52 0.120+^ 

1.92-1.96 0.220lg;jy 2.16-2.20 0.408±g;^| 2.52-2.68 0.050tg;g^ 

1.96-2.00 0.136lo;J?5 2.20-2.24 0.070^;^ 2.68-2.84 0.026tg;^ 

2.00-2.04 0.446^;^ 2.24-2.28 0.1741^ 2.84-3.00 0.026±g^ 

2.04-2.08 0.315lo;n3 2.28-2.36 0.069^:030 3.00-3.16 0.032tg;g^ 

2.08-2.12 0.519lo;l!? 2.36-2.44 0.1121°'°^ 3.16-3.32 0.012±°°^ 




FIG. 20: Cross section for the reaction e+e"^0(lO2O)/ o (98O) 
as a function of effective e + e~ cm. energy obtained from the 
Jf 4 



drop near 2.2 GeV/c 2 , but the statistical errors are large 
and no conclusion can be drawn from this mode alone. 
Possible interpretations are discussed in Section VII. 




X (4K) 



FIG. 21: Distribution of \ 2 from the three-constraint fit for 
K + K~ K + K~ candidates in the data (points). The open 
histogram is the distribution for simulated signal events, nor- 
malized as described in the text. The hatched histogram rep- 
resents the background from non-ISR events, estimated as 
described in the text. The cross-hatched histograms is for 
simulated ISR K + K~tt + tt~ events. 



VI. THE K+K~ K+K~ FINAL STATE 
A. Final Selection and Background 

Figure 21 shows the distribution of x\k for the 
K + K~K + K~ candidates as points, and the open his- 
togram is the distribution for simulated K + K~K + K~ 
events, normalized to the data in the region x\k < 5 
where the backgrounds and radiative corrections are 
small. The hatched histogram represents the back- 
ground from e + e~ — > qq events, evaluated as for the 
other modes. The cross-hatched histogram represents 
the background from simulated ISR K + K ~ir + ir~ events 
with both charged pions misidentified as kaons. 

We define signal and control regions of x\k < 20 and 
20 < x\k < 40, respectively. The signal region contains 
2,305 data and 20,616 simulated events, and the con- 
trol region contains 463 data and 1,601 simulated events. 



Figure 22 shows the K + K~ K + K~ invariant mass distri- 
bution from threshold up to 5 GeV/c 2 for events in the 
signal region as points with errors. The qq background 
(hatched histogram) is small at low masses, but dominant 
above about 4.5 GeV/c 2 . Since the ISR K + K~ Tr + Tr~ 
background does not peak at low x\k values, we in- 
clude it in the background evaluated from the control 
region, according to the method explained in Sec. IV A. 
It dominates this background, which is 10% or lower at 
all masses. The total background is shown as the open 
histogram in Fig. 22. 

We subtract the sum of backgrounds from the number 
of selected events in each mass bin to obtain a number 
of signal events. Considering uncertainties in the cross 
sections for the background processes, the normalization 
of events in the control region, and the simulation statis- 
tics, we estimate a systematic uncertainty on the signal 
yield of less than 5% in the 2-3 GeV/ c 2 region, increasing 
to 10% in the region above 3 GeV/c 2 . 
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4 5 
m(K + K"K + K") (GeV/c 2 ) 

FIG. 22: Invariant mass distribution for K + K~ K + K~ candi- 
dates in the data (points). The hatched and open histograms 
represent, cumulatively, the non-ISR background and the ISR 
background from the control region, which is dominated by 
the contribution from ISR K + K~ tt + -k~ events. 




m(K K'K K") (GeV/c ) 



FIG. 23: (a) Invariant mass distributions for simulated 
K + K~ K + K~ events in the signal (open) and control 
(hatched) regions (see Fig. 21); (b) net reconstruction and 
selection efficiency as a function of mass obtained from this 
simulation (the curve represents a 3 rd order polynomial fit). 



interval by the number of generated ones in that interval 
to obtain the efficiency shown as the points in Fig. 23(b). 
It is quite uniform, and we fit a third order polynomial, 
which we use to extract the cross section. A factor of 0.7 
is again applicable for only hadronic system efficiency due 
to the limited angular coverage of the ISR photon simula- 
tion. A simulation assuming dominance of the (f)K + K~ 
channel, with the K + K~ pair in an S-wave, gives con- 
sistent results, and we apply the same 5% systematic 
uncertainty as for the other final states. We correct for 
mis-modeling of the track finding and kaon identifica- 
tion efficiencies as in Sec. IV B, and for the shape of the 
xIk distribution analogously, using the result in Ref. [11], 
(3.0 ±2.0)%. 



K+K K+K 

We calculate the e+e~^> K+K~K+K~ cross section 
in 40 MeV £ c . m . intervals from the analog of Eq. 2, us- 
ing the invariant mass of the K+K~K+K~ system to 
determine the effective cm. energy. We show this cross 
section in Fig. 24 and list it in Table IX. It rises to a 
peak value near 0.1 nb in the 2.3-2.7 GeV region, then 
decreases slowly with increasing energy. The only sta- 
tistically significant narrow structure is the J/ip peak. 
Again, dC includes corrections for vacuum polarization 
that should be omitted from calculations of — This 
supersedes our previous result [11]. 




B. Selection Efficiency 



(GeV) 



FIG. 24: The e+ e~ ^ K+ K~ K+ K~ cross section as a func- 
tion of the effective e + e~ cm. energy measured with ISR data 
at BABAR. The errors are statistical only. 



The detection efficiency is determined as for the other 
two final states. Figure 23(a) shows the simulated 
K+K~ K+K~ invariant-mass distributions in the signal 
and control regions from the phase space model. We di- 
vide the number of reconstructed events in each mass 



The simulated K+K~K+K~ invariant mass resolu- 
tion is 3.0 McV/c 2 in the 2.0-2.5 GeV/c 2 range, increas- 
ing with mass to 4.7 MeV/c 2 in the 2.5-3.5 GeV/c 2 range. 
Since the cross section has no sharp structure except for 
the J/tp peak, we again make no correction for resolution. 
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TABLE IX: Measurements of the e + e — * K + K K + K cross section (errors are statistical only). 



Ec m . (GeV) 


a 


(nb) 


Ecm. (GeV) 


a 


(nb) 


Ec.m. (GeV) 


a 


(nb) 


E c . m . (GeV) 


a (nb) 


2.02 


0.006 


± 


0.004 


2.66 


0.075 


± 


0.015 


3.30 


0.025 


± 


0.009 


3.94 


0.012 ± 0.006 


2.06 


0.011 


± 


0.006 


2.70 


0.102 


± 


0.015 


3.34 


0.025 


± 


0.009 


3.98 


0.012 ± 0.006 


2.10 


0.019 


± 


0.007 


2.74 


0.069 


± 


0.014 


3.38 


0.024 


± 


0.010 


4.02 


0.010 ± 0.006 


2.14 


0.034 


± 


0.011 


2.78 


0.063 


± 


0.014 


3.42 


0.034 


± 


0.010 


4.06 


0.009 ± 0.005 


2.18 


0.040 


± 


0.013 


2.82 


0.051 


± 


0.012 


3.46 


0.036 


± 


0.009 


4.10 


0.006 ± 0.005 


2.22 


0.087 


± 


0.016 


2.86 


0.024 


± 


0.011 


3.50 


0.032 


± 


0.009 


4.14 


0.008 ± 0.006 


2.26 


0.064 


± 


0.018 


2.90 


0.054 


± 


0.012 


3.54 


0.025 


± 


0.009 


4.18 


0.011 ± 0.005 


2.30 


0.082 


± 


0.017 


2.94 


0.045 


± 


0.013 


3.58 


0.031 


± 


0.009 


4.22 


0.011 ± 0.005 


2.34 


0.079 


± 


0.018 


2.98 


0.045 


± 


0.011 


3.62 


0.014 


± 


0.007 


4.26 


0.012 ± 0.005 


2.38 


0.084 


± 


0.017 


3.02 


0.063 


± 


0.013 


3.66 


0.019 


± 


0.008 


4.30 


0.004 ± 0.005 


2.42 


0.070 


± 


0.016 


3.06 


0.049 


± 


0.012 


3.70 


0.028 


± 


0.008 


4.34 


0.003 ± 0.006 


2.46 


0.092 


± 


0.018 


3.10 


0.287 


± 


0.024 


3.74 


0.008 


± 


0.005 


4.38 


0.013 ± 0.005 


2.50 


0.082 


± 


0.015 


3.14 


0.050 


± 


0.012 


3.78 


0.026 


± 


0.007 


4.42 


0.012 ± 0.005 


2.54 


0.091 


± 


0.015 


3.18 


0.054 


± 


0.011 


3.82 


0.020 


± 


0.007 


4.46 


0.006 ±-0.004 


2.58 


0.077 


± 


0.017 


3.22 


0.049 


± 


0.011 


3.86 


0.007 


± 


0.006 


4.50 


-0.001 ±-0.004 


2.62 


0.095 


± 


0.015 


3.26 


0.045 


± 


0.010 


3.90 


0.012 


± 


0.006 


4.54 


0.000 ± 0.004 



The errors shown in Fig. 24 and Table IX are statisti- 
cal only. The point-to-point systematic errors are much 
smaller than this, and the errors on the normalization are 
summarized in Table X, along with the corrections ap- 
plied to the measurement. The total correction is ±10%, 
and the total systematic uncertainty is 9% at low mass, 
increasing to 13% above 3 GeV/c 2 . 



TABLE X: Summary of corrections and systematic uncer- 
tainties on the e + e~— > K + K~ K + K~ cross section. The total 
correction is the linear sum of the components and the total 
uncertainty is the sum in quadrature. 



Source 


Correction 




Uncertainty 


Rad. Corrections 




1% 




Backgrounds 




5% 


m iK < 3 GeV/c 2 






10% 


m 4 j f > 3 GeV/c 2 


Model Dependence 




5% 




X4K Distribution 


±3% 


2% 




Tracking Efficiency 


±3% 


2% 




Kaon ID Efficiency 


±4% 


4% 




ISR Luminosity 




3% 





Total 



±10% 



9% 
13% 



m4K < 3 GeV/c 
iri4K > 3 GeV/c 



D. The 4>(1020)K + K " Intermediate State 

Figure 25 shows the invariant mass distribution for all 
K + K~ pairs in the selected K + K~ K + K~~ events (4 en- 
tries per event) as the open histogram. A prominent 
4> peak is visible, along with possible peaks at 1.5, 1.7 
and 2.0 GeV/c 2 . The hatched histogram is for the pair 
in each event with mass closest to the nominal <f> mass, 
and indicates that the <pK + K~ channel dominates the 
K + K~ K + K~ final state. Our previous finding of very 
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FIG. 25: Invariant mass distributions for all K + K~ pairs in 
selected e + e~— > K + K~ K + K~ events (open histogram) and 
for the combination in each event closest to the </>-meson mass 
(hatched). 



little <j) signal [11] was incorrect due to an error in the 
analysis algorithm. 

If the pair mass closest to the <fi mass is within 
10 MeV/c 2 of the </> mass, then we include the invariant 
mass of the other K + K~ combination in Fig. 26. The 
contribution from events in the J/ ip peak is shown as the 
hatched histogram which is in agreement with the BES 
experiment [24] which studied the structures around 1.5, 
1.7 and 2.0 GeV/c 2 in detail. There is no evidence for 
the 4>fo channel, but there is an enhancement at thresh- 
old that can be interpreted as the tail of the /o(980). 
This is expected in light of the 0/q cross sections mea- 
sured above in the K + K~ir + n~ 
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1.8 2 2.2 2,4 
m(K + K") (GeV/c 2 ) 

FIG. 26: Invariant mass distribution for K + K~ pairs in events 
in which the other K + K~ pair has mass closest to and within 
10 MeV/c 2 of the nominal <f> mass (open histogram). The 
hatched histogram is for the subset with a K + K~ K + K~ 
mass in the J/tp peak. 




2.8 3 
E c . m . (GeV) 

FIG. 27: The e+e"^ <j(>(1020)/o(980) cross section measured 
in the K + K~-k + -k~ (circles) and K + K~ tt°it (squares) final 
states. The hatched histogram shows the simulated cross sec- 
tion, assuming no resonant structure. The solid (dashed) line 
represents the result of the one-resonance (no-resonance) fit 
described in the text. 



states. However the statistics and uncertainties in the 
/o(980) — » K + K~ lincshape do not allow a meaningful 
extraction of the cross section in this final state. 

We observe no significant structure in the K + K^K ± 
mass distribution. We use these events to study the pos- 
sibility that part of our <fiir + ir~ signal is due to tfiK + K~ 
events with the two kaons not from the <p taken as pions. 
No structure is present in the reconstructed K + K~tt + it~ 
invariant mass distribution from these events. 



VII. 



<t>f NEAR THRESHOLD 



The behavior of the e + e _ — > 4>fo cross section near 
threshold shows a structure near 2150 MeV/c 2 , and we 
have published this result in Ref. [14]. Here we pro- 
vide a more detailed study of this cross section in the 
1.8-3 GeV region. In Fig. 27 we superimpose the cross 
sections measured in the K + K~iy + tt~ - and K+K-n°n 
final states (shown in Figs. 13 and 20); they are con- 
sistent with each other. The K + K~ K + K~ cross section 
(Sec. VI D) is also consistent with the presence of a struc- 
ture near 2150 MeV/c 2 and shows a contribution from 
the (f>fo channel, but since we cannot extract a meaning- 
ful </>/o cross section, we do not discuss this final state 
further. 

First, we attempt to reproduce this spectrum with a 
smooth threshold function. In the absence of resonances, 
the only theoretical constraint on the cross section well 
above threshold is that it should decrease smoothly with 



increasing 



E, 



However the form of the cutoff at 



diate resonances and the final state particle spins, phase 
space and detector resolution. The model discussed in 
Sec. IV F takes the <f) and /o(980) lineshapes, the spins 
of all particles and their phase space into account, and 
postulates a simple E~^ dependence of the cross section. 
For the e + e~^ <f)fo reaction, it predicts the cross section 
shown as the hatched histogram in Fig. 27, normalized 
to the same total area as the data. It shows a sharp rise 
from the threshold with a peak near 2070 MeV and is 
inconsistent with the data. 

To account for uncertainties in the fo width and the 
shape of the cross section well above threshold, we seek a 
functional form that describes the simulation and whose 
parameters can be varied to cover a reasonable range of 
possibilities. This can be achieved by the product of a 
phase space term, an exponential rise and a second order 
polynomial: 

ow(m) = -P(m) ' A nr (n), (4) 
A nr {ij) = ao ■ (1 - e~^' a ^) • (1 + a 2l i + a zl x 2 ), 



PQi) = y / l-m 2 1 /(m + /i) 2 I M 



E r 



m , 



threshold is determined by the properties of the interme- 



where the ai are free parameters, ao is a normalization 
factor, and P(/x) is a good approximation of the two-body 
phase space for particles with similar masses. Both the 
0(1020) and /o(980) have small but finite widths, and 
our selection criterion of m(7T7r) > 0.85 GeV/c 2 defines 
an effective minimum mass, mo = 1.8 GeV/c 2 . Fitting 
Eq. 4 to the simulated cross section yields the a% values 
listed in the first column of Table XI. Fitting to the data 
with all on fixed to these values and ctq floating yields 
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X 2 /n.d.f.= 86/(56 — 2). Floating a-i and 03 in addition, 
we obtain % 2 /n.d.f.= 85/(56 — 4) with a confidence level 
(C.L.) of 0.0025. If we float all parameters in Eq. 4, 
the fit yields x 2 /n.d.f.= 80/(56- 5) with C.L. of 0.0053. 
The results of these fits are listed in Table XI, and the 
latter is shown as the dashed curve on Fig. 27; all fits are 
inconsistent with the data. 



TABLE XI: Parameter vaiues, x 2 values and confidence levels 
from the fits of Eq. 4 to the data described in the text. An 
asterisk denotes a value that was fixed in that fit. 



Fit 


All a, fixed 


Only 01 fixed 


All ai free 


ffo 
ai 

a 2 

a-i 

X 2 /n.d.f. 

P(x 2 ) 


1.19±0.03 
0.218* 
-1.68* 
0.81* 
86.4/54 
0.0035 


1.23±0.03 
0.218* 
-1.51±0.15 

0.66±0.14 
85.3/52 
0.0025 


1.09±0.01 
0.174±0.012 
-1.49±0.12 
0.63±0.11 
80.5/51 
0.0053 



We now add a resonance and fit the data with the 
function 

aiM = • \A n Me^ + A rl (p)\ 2 , (5) 

m l ~ ^c.m. — ^c.m.l 1 

where mi and Ti are the mass and width of the reso- 
nance, ai is its peak cross section, and ipi is its phase rela- 
tive to the non-resonant component. We obtain good fits 
both assuming no interference between the two compo- 
nents, ipi =7r, and with tpi floating. The result of the lat- 
ter fit is shown as the solid curve on Fig. 27. The data are 
somewhat above this curve near 2.4 GeV/c 2 and a fit with 
two resonances can also describe the data. Due to the 
sharp drop near 2.2 GeV/c 2 , the single-resonance fit with 
interference gives a resonance mass about 30 McV/c 2 
higher than the other two fits. All these fits, with or with- 
out resonances, give a peak non-resonant cross section in 
the range 0.3-0.4 nb, which is of independent theoretical 
interest, because it can be related to the (f> — > /o (980)7 
decay studied at the 0- factory [25]. 

Under the hypothesis of one resonance interfering with 
the non-resonant component, the fit gives the resonance 
parameters 

01 = 0.13 ±0.04 nb, 
mi = 2.175 ±0.010 GeV/c 2 , 
Ti = 0.058 ±0.016 GeV, and 
■01 = -0.57 ± 0.30 radians, 

along with x 2 /n.d.f.= 37.6/(56- 9) (C.L. 0.84). We can 
estimate the product of its electronic width and branch- 
ing fraction to 0/o as 

' T ee - = (2.5 ± 0.8 ± 0.4) eV , 



where wc fit the product TiUi to reduce correlations, and 
the conversion constant C — 0.389 mb(GeV/c 2 ) 2 . The 
second error is systematic and corresponds to the nor- 
malization errors on the cross section. 

The significance of the structure calculated from the 
change in x 2 between the best fit and the null hypothe- 
sis is 6.2 standard deviations. Since this calculation can 
be unreliable in the case of low statistics and functions 
that vary rapidly on the scale of the bin size, we per- 
form a set of simulations in which we generate a number 
of events according to a Poisson distribution about the 
number observed in the data and with a mass distribu- 
tion given by either the simulation or fitted function in 
Fig. 27 without resonant structure. On each sample, we 
perform fits to Eqs. 4 and 5 and calculate the difference 
in x 2 - The fraction of trials giving a x 2 difference larger 
than that seen in the data corresponds to a significance 
of approximately 5 standard deviations. 

We search for this structure in other submodes with 
different and/or fewer intermediate resonances. The to- 
tal cross sections are dominated by K*Kn channels, and 
the K *°K + 7r~ cross section is shown in Fig. 8. There 
is no significant structure in the 2.1-2.5 GeV region, but 
the point-to-point statistical uncertainties are large. If wc 
remove events within the bands in Figs. 5 and 18, then 
most of the events containing a K* are eliminated and 
we obtain the raw mass distributions shown as the points 
with errors in Figs. 28 and 29, respectively. Both distri- 
bution show evidence of a structure around 2.15 GeV/c 2 
and the K + K ~tt + tt~ distribution also shows a structure 
near 2.4 GeV/c 2 . We cannot exclude the presence of these 
structures in events with a K* , but we can conclude that 
they do not dominate those events, whereas they com- 
prise a substantial fraction of the remaining events in 
that mass region. 

Applying the further requirement that the dipion mass 
be in the range 0.85-1.10 GeV/c 2 , we remove most of 
the events without an f , and obtain the mass distribu- 
tions shown as the hatched histograms in Figs. 28 and 29. 
Peaks are visible at both 2.15 GeV/c 2 and 2.4 GeV/c 2 in 
both distributions, and they contain enough events to 
account for the corresponding structures in the distribu- 
tions for all non-K* events. These peaks contain at least 
as many events as are present in the (fifo samples, but 
the non-resonant components are higher and there is a 
substantial kinematic overlap between K + K~f events 
and K*Ktt events in this mass range. 

Since this /o(980) band appears to contain a large frac- 
tion of the events within the structure, we now consider 
all selected events with a dipion mass inside or outside 
this range. Figure 31 shows the mass distribution for all 
selected K + K~it tt events as the open histogram, and 
the subsets of events with 7r°7r° mass inside and outside 
the range 0.85-1.10 GeV/c 2 as the hatched and cross- 
hatched histograms, respectively. It is evident that the 
K + K~f channel contains the majority of the structure 
in the 2.0-2.6 GeV/c 2 range. 

We show the corresponding distributions for the 
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FIG. 28: Invariant mass distribution for all selected 
K + K~ tt + tt~ events lying outside the _ft"" (892) bands of 
Fig. 5 (points), and the subset of these events with 0.85 < 
m{TT + 7v-) < 1.10 GeV/c 2 (hatched). 
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FIG. 29: Invariant mass distribution for all selected 
K + K'TT°n° events lying outside the 7f*°(892) bands of 
Fig. 18 (points), and the subset of these events with 0.85 < 
m(7rV) < 1.10 GeV/c 2 (hatched). 




3.5 4 4.5 5 
m(K + K"7tV) (GeV/c 2 ) 

FIG. 30: The K + K ~tt + tt~ invariant mass distribution for all 
selected events (open histogram), and for those with a tv + tt~ 
mass inside (cross-hatched) or outside (hatched) the fo band 
as defined in the text. 




3.5 4 4.5 5 
m(K + K"3lV) (GeV/c 2 ) 

FIG. 31: The K + K~ 7r°7r° invariant mass distribution for all 
selected events (open histogram), and for those with a 7r°7r° 
mass inside (cross-hatched) or outside (hatched) the fo band 
as defined in the text. 



K + K~iy + tt~ events in Fig. 30. Due to the presence of 
the p°, the relative fo contribution is much smaller in this 
final state, but the events in the fo band show clear in- 
dications of structure in the 2.0-2.4 GeV/c 2 region. The 
remaining events may also have structure in this region, 
but the statistical significance is marginal and it could 
be due to other sources, such as the </>/ 2 (1270) threshold 
at 2.3 GeV/c 2 . 

Figures 32 and 33 show enlarged views of the mass 
distributions within the f bands from Figs. 30 and 31, 
respectively. The two-peak structure is more evident 



here than in the <j>fo events. The 0.85 < m(7T7r) < 
1.10 GeV/c 2 requirement gives enough phase space for 
K + K~ invariant mass to cover the region from threshold 
to -1.3 GeV/c 2 for m(K+K-nn) w 2.15 GeV/c 2 . From 
the measured kaon form factor we expect to find only 
about two-thirds of K + K~ P-wave in our fitted <j> peak. 
Since the non-ISR and ISR 7T7r7r7r backgrounds have not 
been subtracted and the samples contain an unknown 
mixture of intermediate states, we fit them with a modi- 
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TABLE XII: Summary of parameters obtained from the fits described in the text to the K + K tt + tt and K + K 7r l V° events 
with dipion mass in the /o(980) band. An asterisk denotes a value that was fixed in that fit. 



No Resonance One Resonance Two Resonances 



Fit 


K+R-TT + IT- 


K+K-it tv 




K+K-n°ir° 


K+K-TT+TY- 




N nr 


7204±775 


991±202 


8466±334 


722±112 


6502±476 


117±89 


ai 


0.f8f±0.012 


0.134±0.017 


0.224±0.024 


0.197±0.048 


0.201±0.035 


0.143±0.053 


a 2 


-0.75±0.21 


-1.47±0.38 


-0.89±0.17 


-0.36±0.10 


-0.44±0.15 


5.80±2.36 


a3 


0.09±0.17 


0.75±0.35 


0.17±0.08 


-0.28±0.14 


-0.15±0.12 


-5.26±1.75 


(24 


0.75* 


0.50* 


0.75* 


0.50* 


0.75* 


0.50* 


Ni 


0* 


0* 


116±95 


149±36 


163±70 


192±44 


mi ( GeV/c 2 ) 






2.192±0.014 


2.169±0.020 


2.187±0.013 


2.154±0.029 


ri ( GeV) 






0.071±0.021 


0.102±0.027 


0.066±0.018 


0.110±0.022 


ipi (rad) 






-0.60±0.41 


-1.02±0.19 


-1.10±0.14 


-1.04±0.23 


N 2 


0* 


0* 


0* 


0* 


16±16 


6±5 


m 2 (GeV/c 2 ) 










2.47±0.07 


2.45±0.04 


r 2 ( GeV) 










0.077±0.065 


0.062±0.102 


ip2 (rad) 










0.28±1.06 


1.41±1.29 


xVn-d.f. 


62.8/41 


38.1/21 


35.6/37 


13.0/17 


31.4/34 


9.7/13 


P(x 2 ) 


0.016 


0.012 


0.54 


0.74 


0.60 


0.72 



fied version of Eq. 5 that allows up to two resonances, 

F(jjt) = (a 4 -A nr ) 2 (6) 
+ \{l-a 4 )A nr + A rl e^ 1 +A r2 e^ 2 \ 2 . 

Here, the normalization is in terms of events rather than 
cross section (<7j — ► Ni) and a fraction a 4 of the non- 
resonant component does not interfere with the reso- 
nances. We first fit the distribution with no resonances 
(and a 4 = 1). The results are shown as the dashed lines 
in Figs. 32 and 33 and listed in Table XII; both are in- 
consistent with the data. 

We next include one resonance in the fit. The param- 
eter a 4 is not well constrained by the data and its value 
has a small influence on all other fit parameters except 
for the number of events assigned to the resonance, so we 
present results with a 4 fixed to the reasonable values of 
0.75 and 0.50 for the K+K'-k+'k- and K+K-tt tt data, 
respectively. The results are shown as the solid lines in 
Figs. 32 and 33 and listed in Table XII. The fit quality is 
good in both cases, the fitted resonance parameters are 
consistent with those from the (j)fo study, and the cal- 
culated significance of the structure for the K + K ~tt + it~ 
data is similar, 5.2 standard deviations. The K + K -7r°7r° 
data show much more pronounced structure than in the 
cj>fo study, allowing a full fit to this sample with a signif- 
icance of 5.0 standard deviations. 

We then add a second resonance to the fit, keeping a 4 
fixed and floating all other parameters. The results are 
shown as the dotted lines in Figs. 32 and 33, and listed 
in Table XII. These fits are also of good quality, but 
do not change the % 2 CL or the parameters of the first 
resonance significantly. We also perform fits with no in- 
terference between the non-resonant component and any 
resonance (04 — 1), obtaining good quality fits for both 
one resonance and two resonances with relative phase 



7r/2. The fitted resonance parameters are consistent in 
all cases, except that the mass of the first resonance is 
lower by about 50 MeV/c 2 , similar to the 30 MeV/c 2 shift 
seen in the <j>fo study. 

From these studies we conclude that we have observed 
a new vector structure at a mass of about 2150 MeV/c 2 
with a significance of over six standard deviations. It 
decays into K + K~ fo(980), with the K + K~ pair pro- 
duced predominantly via the </>(1020). There is an ad- 
ditional structure at about 2400 MeV/c 2 , and the two 
structures can be described by either two resonances or 
a single resonance that interferes with the non-resonant 
K + K _ Jo (980) process. More data and searches in other 
final states are needed to understand the nature of these 
structures. 

If the main structure is due to a resonance, then it 
is relatively narrow and might be interpreted as the 
strange analog of the recently observed charmed Y(4260) 
state [6], which decays to J/tpir + ir~. The value of 
B<f,f ■ r ee = (2.5 ± 0.8 ± 0.4) eV measured here is similar 
to the value of B Y ^j/^+ n - ■ T^ e = (5.5 ± 1.0 ± 0.8) eV 
reported in Ref. [6]. 



VIII. THE CHARMONIUM REGION 

The data at masses above 3 GeV/c 2 can be used to 
measure or set limits for the branching fractions of nar- 
row resonances, such as charmonia, and the narrow J/ip 
and ip{2S) peaks allow measurements of our mass scale 
and resolution. Figures 34, 35 and 36 show the in- 
variant mass distributions for the selected K + K~n + n~ , 
K + K ~ir°ir° and K + K~K + K~ events, respectively, in 
this region, with finer binning than in the corresponding 
Figs. 2, 15 and 22. We do not subtract any background 
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2.4 2.6 2.8 

m(K + K" f„) (GeV/c 2 ) 

FIG. 32: The K + K~ tv + iv~ invariant mass distribution in the 
K+ K~f b(980) threshold region for events with a 7r + 7r mass 
inside the /o band. The lines represent the results of the 
fits including no (dashed), one (solid) and two (dotted) reso- 
nances described in the text. 




2.4 2.6 2.8 
m(K + K" f ) (GeV/c 2 ) 

FIG. 33: The K + K~ ty°-k° invariant mass distribution in the 
K + K~ /o(980) threshold region for events with a 7r°7r° mass 
inside the /o band. The lines represent the results of the 
fits including no (dashed), one (solid) and two (dotted) reso- 
nances described in the text. 



from the K + K~ir + ir~ or K + K~K + K~ data, since it is 
small and nearly uniformly distributed, but we use the 
XkkttOk con t r °l region to subtract part of the ISR back- 
ground from the K + K~ir°ir° data. Signals from the J/tp 
are visible in all three distributions, and the tp{2S) is vis- 
ible in the K + K~ir + ir~ mode. 

We fit each of these distributions using a sum of two 
Gaussian functions to describe the J/tp and tp(2S) signals 
plus a polynomial to describe the remainder of the dis- 




3.5 3.75 
m(K + K"7tV) (GeV/c 2 ) 

FIG. 34: Raw invariant mass distribution for all selected 
e + e~ ^ K + K~-k + tt~ events in the charmonium region. The 
line represents the result of the fit described in the text. 



tribution. We take the signal function parameters from 
the simulation, but let the overall mean and width float 
in the fit, along with the amplitude and the coefficients 
of the polynomial. The fits are of good quality and are 
shown as the curves on Figs. 34, 35 and 36. In all cases, 
the fitted mean value is within I MeV/c 2 of the PDG [5] 
J '/tp or ip(2S) mass, and the width is consistent within 
10% with the simulated resolution discussed in Sec. IV C, 
VC or VIC. 

The fits yield 1586 ± 58 events in the J/tp peak for the 
K+K-jt+tt- final state, 203 ± 16 events for K+R-^^ 
and 156± 15 events for K + K~ K + K~ . From these num- 
bers of observed events in each final state /, Nj/^f, we 
calculate the product of the J/tp branching fraction to / 
and the J/tp electronic width: 



B 



J/iHf ' 1 ee 



■£ J H> 



N J/4>->f ■ m2 



J/-4> 



67t 2 • dC/dE ■ tf(mj/^) ■ C 



(J) 



where dC/dE = 89.8 nb~ 1 /MeV and ef(mj/^) are the 
ISR luminosity and corrected selection efficiency, respec- 
tively, at the J/tp mass and C is the conversion constant. 
We estimate e K + K - 7r + 7r - =0.202, e K + K - n o n o =0.069 and 
ck+k-k+k- =0.176. 

Using Til 4 ' = 5.40 ± 0.18 keV [5], we obtain the 
branching fractions listed in Table XIII, along with 
the measured products and the current PDG val- 
ues. The systematic errors include a 3% uncertainty 
on Tie . The branching fractions to K + K~ir + ir~ 
and K + K~ K + K~ arc more precise than the current 
PDG values, which were dominated by our previous 
results of (6.25±0.80)xl0" 3 and (7.4±1.8) x 10~ 4 , re- 
spectively [11]. This is the first measurement of the 
K + K - 7T 7r° branching fraction. 
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m(K + K";i ji ) (GeV/c ) 

FIG. 35: Invariant mass distribution for e + e - ^ K + K~ ir ^ 
events in the charmonium region, after partial background 
subtraction. The line represents the result of the fit described 
in the text. 
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FIG. 36: Raw invariant mass distribution for all selected 
e + e~ —> K + K~ K + K~ events in the charmonium region. The 
line represents the result of the fit described in the text. 



These fits also yield 91±15 K + K~ir + ir~ events in the 
ip{2S) peak, but no other significant signals. We expect 
6.3 events from ip(2S) — > J / ipir + n~ K + K~ir + ir~ from 
the relevant branching fractions [5], which is less than 
the statistical error. Subtracting this contribution and 
using a calculation analogous to Eq. 7, with dC/dE = 
115.3 nb _1 / MeV ! we obtain the product of the ip(2S)^> 
K + K~7t + Tr~ branching fraction and its electronic width. 



Dividing by the world average value of T 



[5], 



As noted in Sec. IV D and shown in Fig. 5, 
the K + K~tt + tt~ final state is dominated by the 
K*°(892)Kir channels, with a small fraction seen in the 
K*°(892)Kf(U30) + c.c. channels. Figure 37 shows a 
scatter plot of the invariant mass of a J fiT ± 7r =F pair ver- 
sus that of the K + K ~7r + 7r~ system in events with the 
other K^ir pair near the K*°(892) mass, i.e. within 
the bands in Fig. 5(a) with overlapped region taken only 
once. There is a large concentration of entries in the 
J/ip band with K^if* masses near 1430 MeV/c 2 , but 
no solid evidence for a horizontal band corresponding 
to i^2°(1430) production other than in J/ip decays. We 
show the X ± 7r T mass projection for the subset of events 
with a K + K~ir + ir~ mass within 50 MeV/c 2 of the J/tp 
mass in Fig. 38 as the open histogram. The hatched his- 
togram is the projection for events with a K + K ~t: + tt~ 
mass between 50 and 100 MeV/c 2 below the J/tp mass. 




obtain the branching fraction listed in Table XIII; it is 
consistent with the current PDG value [5]. 



m(K + K"7tV) (GeV/c 2 ) 

FIG. 37: The K if^ invariant mass versus K + K~ tt + tt~ in- 
variant mass for events with the other K^-k^ combination in 
the A"" (892) bands of Fig. 5(a). The overlapped region is 
taken only once. 

The J/ip component appears to be dominated by the 
Kf(U30). Also seen is_a small signal from K*°(892) 
indicating the K*°(892)K*°(892) decay of J/ip: this is 
also seen as an enhancement in the vertical J/ip band in 
Fig. 37. The enhancement at 1.8 GeV/c 2 of Fig. 38 can be 
explained by the J/ip decay into K*°(892)K 2 {1770) + c.c. 
(or K*° (892) K 2 {1820) + c.c.) , a mode which has not pre- 
viously been reported. Subtracting the number of side- 
band events from the number in the J/ip mass window, 
we obtain 317±23 events with a K it^ mass in the range 
1200-1700 MeV/c 2 , which we take as a measure of J/ip 
decays into K*° (892)K* 2 ° (WiQ) , 25 ± 8 events in the 
0.8-1.0 GeV/c 2 window for the K*°(892)K*°(892) de- 
cay and 110 ± 14 events for the X*°(892)X 2 (1770) or 
K* (892) J ftr 2 (1830) final state in the 1.7-2.0 GeV/c 2 re- 
gion. We convert these to branching fractions using Eq. 7 
and dividing by the known branching fractions of excited 
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FIG. 38: The K^tt^ mass projection from Fig. 37 for events 
with a K + K~ty + tt~ mass within 50 MeV/c 2 of the J/tp mass 
(open histogram) and 50-400 MeV/c 2 below (hatched). 



kaons [5]. The results are listed in Table XIII: they 
are considerably more precise than the PDG values. We 
cannot calculate Bj/^_ >k*°k 2 (i770) because no branching 
fractions of i^2(1770) or ^(1830) to Kit are reported. 

We study decays into <pir + ir~ and 07r°7r° using the 
mass distributions shown in Figs. 39 and 40, respectively. 
The open histograms are for the events with a K + K~ 
mass within the (p bands of Figs. 9(c) and 19(c). The 
cross-hatched histogram in Fig. 39 is from the </> side- 
bands of Fig. 9(c) and represents the dominant back- 
ground in the </>7r + 7r~ mode. The hatched histogram in 
Fig. 40 is from the Xkk^Ott con trol region and represents 
the dominant background in the 07r°7r° mode. Subtract- 
ing these backgrounds, we find 103±12 J/tp — ► <pir + ir~ 
events, 23±6 J/tp -> (pn n events, and 10±4 tp(2S) -> 
(f>ir + Tr~ events. We convert these to branching fractions 
and list them in Table XIII. This is the first measure- 
ment of the J/tp — > 07r°7r° branching fraction, and the 
other two are consistent with current PDG values. 

If the Y (4260) has a substantial branching fraction into 



</>7r + 7r~, then we would expect to see a signal in Fig. 39. 
In the mass range |m(</>7r + 7r~) — m(Y)\ < 0.1 GeV/c 2 , 
we find 10 events, and assuming a uniform distribu- 
tion we estimate 9.2 background events from the 3.8- 
5.0 GeV/c 2 region. This corresponds to a signal of 
0.8 ± 3.3 events or a limit of < 5 events at the 90% C.L. 
Using d£/dE = 147.7 nb' 1 / MeV at the F(4260) mass, 
we calculate 8^^+^- • Fj e < 0.4 eV which is well below 
the value of B Y ^j/^ +7r - ■ T^ e = (5.5 ± 1.0 ±0.8) eV [6]. 
No y(4260) signal is seen in any other mode studied here. 

Figures 41(a) and 42 show the corresponding mass dis- 
tributions for (j)fo events, i.e. the subsets of the events in 
Figs. 39 and 40 with a dipion mass in the range 0.85- 
1.10 GeV/c 2 . Signals at the J/tp mass are visible in 
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39: Raw invariant mass distributions for candidate 
e + e~ — > (f>iT + TY~ events (open histogram) and events in the 
<j> side bands of Fig. 9(c) (cross-hatched) in the charmonium 
region. The vertical lines indicate the region used for the 
^(4260) search. 
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40: Raw invariant mass distributions for candidate 
— » (j>-K°-K events (open histogram) and events in the 



Xkktt tt° control region (hatched) in the charmonium region. 



both cases, but <pfo is not the dominant mode of the 
J/tp — > 07r + 7r _ decay. Figure 41(b) shows the tt + tt~ in- 
variant mass distribution for events in the J/tp peak of 
Fig. 39, 3.05 < m(K+K -tt+ tt" )< 3.15 GeV/c 2 . A two- 
peak structure is visible that can interpreted as due to 
the /o(980) and /2(1270) resonances. Fitting the dis- 
tribution in Fig. 41(b) with a sum of two Breit-Wigner 
functions with parameters fixed to PDG values [5], we 
find 19.5 ± 4.5 J/tp -> <pf events and 44 ± 7 J/tp -> <pf 2 
events. From Fig. 42 we estimate 7.0 ± 2.8 4>fo events in 
the 7r°7r° mode. 
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TABLE XIII: Summary of the J/xj) and tp(2S) branching fractions measured in this article. 



Measured 


Measured 


J/ib or i>(2S) Branching Fraction (10 -3 ) 


Quantity 


Value (eV) 


Calculated, this work 


PDG2006 


r J/lp K! 

Lee *-V/i/>— >iT+K- tt+ti— 


36.3 ±1.3 ±2.1 


6.72±0.24±0.40 


6.2 ±0.7 


pJ/V KJ 

J- ee •Dj/^->K+ K~ tt°-k° 


13.6 ±1.1 ±1.3 


2.52±0.20±0.25 


no entry 


-pJ/V r> 

J- ee ■Dj/$-*K+ K~ K+ K~ 


4.11±0.39±0.30 


0.76±0.07±0.06 


0.78 ±0.14 




7.3 ±0.5 ±0.6 


2.7 ±0.2 ±0.2 


6.7 ±2.6 


Tee ■Bj/^^K*°K*Q ' ^K*°^Ktt ' ^KtO^Kir 


0.57±0.18±0.05 


0.11±0.04±0.01 


<0. 5 at 90% C.L. 


■pJ/i> D K? K?„ 


o k _|_n *3 -4-0 9 




no entry 


Tee 'Sj/^^07r + 7T- ' K~ 


2.61±0.30±0.18 


0.98±0.11±0.07 


0.94 ±0.15 


Tee -Bj/^^^O^O • B ( f > _ lK + K - 


1.54±0.40±0.16 


0.58±0.15±0.06 


no entry 




0.50±0.11±0.04 


0.28±0.07±0.02 


0.32 ±0.09 (s=1.9) 


vil^-Bj/^^f - B^ K + K - ■B fo ^ 7I o 7I o 


0.47±0.19±0.05 


0.54±0.21±0.05 


0.32 ±0.09 (s=1.9) 


ri£^'Bj/ii>-Kfrf 3 ■B (j> _ K j rK - ■ B f2 ^ n+7: - 


1.12±0.18±0.09 


0.50±0.08±0.04 


<0. 37 at 90% C.L. 


r V>(2S) „ 

-Lee % >~ > i>(2S)^K+ K~ tt+tt- 


2.56±0.42±0.16 


1.2 ±0.2 ±0.08 


0.72 ±0.05 


pV>(2S) R R 

-Lee D V>(2S)-»07r+ir- ' °4>^K+K- 


0.28±0.11±0.02 


0.27±0.11±0.02 


0.113 ±0.029 


Tie 2S ' > ■B^,(2S)^<t>f '®<t>^K+K- ' Bf ->w+n- 


0.17±0.08±0.02 


0.26±0.12±0.03 


0.090 ±0.033 
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FIG. 41: (a) Raw invariant mass distribution for candidate 
</>/o , /o— >7r + 7r _ events (open histogram) and events in the </> 
side bands (cross-hatched) in the charmonium region; (b) the 
tt + tv~ invariant mass distribution for 07r + 7r~ events from the 
J/ip peak of Fig. 39. The line represents the result of the fit 
described in the text. 



Using Eq. 7 and dividing by the appropriate branching 
fractions, we obtain the J/tp branching fractions listed in 



Table XIII. The measurements of B 



in the it + tt 



and 7r u 7r u decay modes of the fo are consistent with each 
other and with the PDG value, and combined they have 
roughly the same precision as listed in the PDG [5] . This 



is the first measurement of B 



and the value is 



> 
o 




consistent with the previous upper limit [5]. We also 



m(<)> f () ) (GeV/c ) 

FIG. 42: Raw invariant mass distribution for candidate <f>fo, 
events (open histogram) and events in the Xkktv°tv 
control region (hatched) in the charmonium region. 



observe 6 ± 3 ip(2S) —*(f>fo> fo~ >n + ir~ events, which we 
convert to the branching fraction listed in Table XIII; it is 
consistent with the PDG value [5], assuming Bf ^ n + n - — 
2/3. 

In the F(4260) region we have 4 events with an esti- 
mated background of about 1 event. This corresponds to 
3±2 events, or a 90% CL upper limit of 5.6 events. We 
convert this to the limits 

BY^^fo ■ ^e e ■ B<t>^K+K- ' %,-►*•+ 71- < 0.14 eV , 

By^fo-T^ < 0.43 eV, 90%C£, 
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which is again much lower than the corresponding quan- 
tity for the y(4260)-> J/tpn+n- decay. 



IX. SUMMARY 

We use the excellent charged particle tracking and 
identification, and photon detection of the BABAR 
detector to fully reconstruct events of the type 
e + e~ — ► je + e~ — > jR + R~ir + ir~ , -fR + R~ir°ir a and 
jR + R~ R + R~ , where the 7 is radiated from the ini- 
tial state e + or e~ . Such events are equivalent to direct 
e+e~ annihilation at an effective cm. energy correspond- 
ing to the mass of the hadronic system, and we study 
the annihilation into these three final states at low E c , m _, 
from their respective production thresholds up to 5 GeV. 
The K + K ~ir + ir~ and R + R~R + R~ measurements are 
consistent with, and supersede, our previous results [11]. 
This is the first measurement of the R + R~w°ir° final 
state, although some of the results were also presented in 
Ref. [14]. 

The systematic uncertainties on the normaliza- 
tion of the e+e- -> R+R'tt+tt-, R+R-jt ?: and 
R+R- R+R- cross sections are 8%, 10% and 9%, re- 
spectively, for E c . m . <3 GeV, and 10%, 14% and 13% in 
the 3-5 GeV range. The obtained cross sections are con- 
siderably more precise than previous measurements and 
cover this low energy range completely, so they provide 
important input to calculations of the hadronic correc- 
tions to the anomalous magnetic moment of the muon 
and the fine structure constant at the Z° mass. 

These final states exhibit complex resonant substruc- 
ture. In the R + R~tt + tt~ mode we measure the cross 
sections for the first time for the specific channels e + e~^ 
R*°(890)R ~7r+, 07r + 7r~ and <f>f . We also observe sig- 
nals for the p°(770), A"i(1270), A"i(1400), Rf{U3Q) 
and /|°(1270) resonances. It is difficult to disentangle 
these contributions to the final state, and we make no 
attempt to do so in this paper. We note that the p° 
signal is consistent with being due entirely to R\ de- 
cays, and the total cross section is dominated by the 
R *°(892)if ~7r + + c.c. channels, but there is no signifi- 
cant signal for e + e-^ir°(892)ir°(892). 

In the R + R~ir°ir a mode we measure cross sections 
for e+e~ — > </>/ and observe signals for the A"* ± (892) 
and R^ (14:30) resonances. Again, the total cross sec- 
tion is dominated by the R*+(892)R+ir° channels, and 
there is no signal for e+e" -> R*+(892)R*-(892). The 
R + R~ir°ir° final state is not accessible to intermediate 
states containing R\ resonances, and we note that the 
cross section is roughly a factor of four smaller over most 
of the range than the R+R-jt+tt- cross section, consis- 
tent with R*Rir dominance with a factor of two isospin 
suppression of the 7r°7r° final state and another factor of 
two for the relative branching fractions of the neutral and 
charged R* to charged kaons. 

In the R+R- R+R- mode wc find e + e~^(f>R + R~ to 
be the dominant channel. With the current data sample 



we can say little about the other R+R- combination, 
except that there is an enhancement near threshold, con- 
sistent with the (j)fo channel, and that in J/ip decays 
there is structure in the 1.5-2.0 GeV region, consistent 
with that observed by BES [24]. 

The 4>fo cross section measured in the R+R --k+tt- 
final state shows structure around 2.15 GeV and pos- 
sibly 2.4 GeV, and the corresponding measurement in 
the R + R ~ir°Tr° final state is consistent, as reported in 
Ref. [14] . Further investigation here reveals consistent re- 
sults in the R+R- R+R- final state and clear signals in 
the R + R~fo channels, with /o^7r + 7r~ and tt°tt°. The 
signals are predominantly from <pfo, but the relaxation 
of the R+R- mass requirement reveals a strong signal 
in the R + R~n°ir final state. This structure can be in- 
terpreted as a strange partner (with c-quarks replaced 
by s-quarks) of the recently observed V(4260), which has 
the analogous decay mode J / ipTT + it~ , or as an ssss state 
that decays predominantly to <frf . 

We also study charmonium decays into these final 
states and their intermediate channels. All nine of the 
J/ip branching fractions and one of the three tp(2S) 
branching fractions listed in Table XIII are as precise or 
more precise than the current world averages. We do not 
observe the V (4260) in any decay mode. In particular, we 
find that the branching fraction for the Y"(4260)^ </>7r + 7r~ 
decay, that a glueball model [8] predicts, is less than one- 
tenth of that to J/ipir + ir~ . 
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